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Apparatus and Method for Providing Robust Error Recovery 
for Errors that occur in a Lossy Transmission Environment 

BACKGROUND OF THE INVENTION 

1. Related Applications 

This application is a continuation of U.S. Patent Application 
Serial No. 09/016,083, filed January 30, 1998, entitled "Source Coding 
to Provide for Robust Error Recovery During Transmission Losses"; 
which is a continuation-in-part of Application Serial No. 09/002,547, 
filed January 2, 1998, entitled "Image-to-Block Mapping to Provide for 
Robust Error Recovery During Transmission Losses", Application 
Serial No. 09/002,470, filed January 2, 1998, entitled "Source Coding to 
Provide for Robust Error Recovery During Transmission Losses"; 
Application Serial No. 09/002,553, filed January 2, 1998, entitled 
"Multiple Block Based Recovery Method to Provide for Robust Error 
Recovery During Transmission Losses"; which are continuations-in- 
part of Application Serial No. 08/956,632, filed October 23, 1997, 
entitled "Image-to-Block Mapping to Provide for Robust Error 
Recovery During Transmission Losses"; Application Serial No. 
08/957,555, filed October 23, 1997 entitled "Source Coding to Provide 
for Robust Error Recovery During Transmission Losses"; and 
Application Serial No. 08/956,870, filed October 23, 1997, entitled 
"Multiple Block Based Recovery Method to Provide for Robust Error 
Recovery During Transmission Losses". Application Serial No. 
09/016,083, filed January 30, 1998, Application Serial No. 09/002,547, 
filed January 2, 1998, Application Serial No. 09/002,470, filed January 2, 
1998, Application Serial No. 09/002,553, filed January 2, 1998, 
Application Serial No. 08/956,632, filed October 23, 1997, Application 
Serial No. 08/957,555, filed October 23, 1997 and Application Serial No. 
08/956,870, filed October 23, 1997 are herein incorporated by reference. 

2. Field of the Invention 

The present invention relates to providing a robust error 
recovery due to data losses incurred during transmission of signals. - 

3. Art Background 
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A number of techniques exist for reconstructing lost data due to 
random errors that occur during signal transmission. However, these 
techniques cannot handle the loss of consecutive packets of data. 
Consecutive loss of packets of data is described in the art as burst error. 
Burst errors result in a reconstructed signal with such a degraded 
quality that it is easily apparent to the end user. Additionally, 
compression methodologies used to facilitate high speed 
communications compound the signal degradation caused by burst 
errors, thus adding to the degradation of the reconstructed signal. An 
example of burst error loss affecting transmitted and/or stored signals 
is seen in high definition television ("HDTV") signals and mobile 
telecommunication applications wherein compression methodologies 
play an important role. 

The advent of HDTV has led to television systems with a much 
higher resolution than the current standards proposed by the National 
Television Systems Committee ("NTSC"). Proposed HDTV signals 
are predominantly digital. Accordingly, when a color television signal 
is converted for digital use it is common that the luminance and 
chrominance signals are digitized using eight bits. Digital 
transmission of color television requires a nominal bit rate of two 
hundred and sixteen megabits per second. The transmission rate is 
greater for HDTV which would nominally require about 1200 
megabits per second. Such high transmission rates are well beyond 
the bandwidths supported by current wireless standards. Accordingly, 
an efficient compression methodology is required. 

Compression methodologies also play an important role in 
mobile telecommunication applications. Typically, packets of data are 
communicated between remote terminals in mobile 
telecommunication applications. The limited number of 
transmission channels in mobile communications requires an 
effective compression methodology prior to the transmission of 
packets. A number of compression techniques are available to 
facilitate high transmission rates. 
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Adaptive Dynamic Range Coding ("ADRC") and the discrete 
cosine transform ("DCT") coding provide image compression 
techniques known in the art. Both techniques take advantage of the 
local correlation within an image to achieve a high compression ratio. 
However, an efficient compression algorithm results in compounded 
error propagation because errors in an encoded signal are more 
prominent when subsequently decoded. This error multiplication 
results in a degraded video image that is readily apparent to the user. 

SUMMARY OF THE INVENTION 

A method for source coding a signal is described. In particular, 
a signal comprising multiple signal elements is processed. Each signal 
element is encoded to form a bitstream. The bits within a given 
bitstream are distributed across different bitstreams. Thus, the 
parameters describing components of the segment elements are 
distributed across the different bitstreams. The distributing steps 
result in error distribution across multiple levels. Therefore, when 
the distributing steps are reversed by the decoder, a burst transmission 
error becomes a distributed set of localized losses. 

Another method is also described for a multiple level shuffling 
process. A signal is defined as multiple levels wherein each level 
comprises a plurality of frames, a plurality of pixels, and a plurality of 
bits. In one embodiment, shuffling occurs on each level and between 
levels. Multiple level shuffling causes burst error loss to be 
distributed across multiple levels thereby facilitating image 
reconstruction of those areas of the image in which the loss occurred. 

BRIEF DESCRIPTION OF THE DRAWINGS 

The objects, features and advantages of the present invention 
will be apparent to one skilled in the art in light of the following 
detailed description in which: 
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Figure 1 generally illustrates the processes of signal encoding, 
transmission, and decoding. 

Figure 2 illustrates one embodiment of a packet structure. 

Figure 3 is a flow diagram illustrating one embodiment of the 
encoding process in accordance with the teachings of the present 
invention. 

Figure 4 is a flow diagram illustrating one embodiment of the 
decoding process in accordance with the teachings of the present 
invention. 

Figure 5 illustrates one embodiment of image-to-block mapping 
in accordance with the teachings of the present invention. 

Figure 5a illustrates one embodiment of a shuffling pattern 
used in image-to-block mapping. 

Figure 6 is an illustration of exemplary complementary and 
interlocking block structures. 

Figures 7a, 7b, 7c, 7d illustrate one embodiment of shuffling 
patterns for Y blocks within a frame set. 

Figure 8 is an illustration of one embodiment of cumulative 
DR distribution for Buffer 0. 

Figure 8a is an illustration of one embodiment of a partial 
buffering process in accordance with the teachings of the present 
invention. 
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Figure 9 illustrates one embodiment of the intra buffer YUV 
block shuffling process in accordance with the teachings of the present 
invention. 

Figure 10 illustrates one embodiment of the intra group VL- 
data shuffling process in accordance with the teachings of the present 
invention. 

Figure 11 illustrates one embodiment of Q code concatenation 
within a 3-block group in accordance with the teachings of the present 
invention. 

Figure 11a illustrates one embodiment of Q code concatenation 
for frame pairs including motion blocks in accordance with the 
teachings of the present invention. 

Figure 12 illustrates one embodiment of pixel data error caused 
by a 1/6 burst error loss. 

Figure 12a illustrates one embodiment of shuffling Q codes and 
distributing Q code bits in accordance with the teachings of the present 
invention. 

Figure 12b illustrates one embodiment of pixel data error caused 
by a 1/6 burst error loss of redistributed Q codes. 

Figure 12c illustrates one embodiment of pixel data error caused 
by a 1 / 6 burst error loss of reassigned Q codes. 

Figure 13 illustrates one embodiment of MIN shuffling in 
accordance with the teachings of the present invention. 

Figure 13a illustrates one embodiment of Motion Flag shuffling 
and of a fixed length data loss in one frame pair. 
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Figure 14 illustrates one embodiment of a modular shuffling. 

Figure 14a illustrates one embodiment of a modular shuffling 
result and the fixed length data loss associated with the modular 
shuffling. 

Figure 14b illustrates an alternate embodiment of a modular 
shuffling result and the fixed length data loss associated with the 
modular shuffling. 

Figure 14c illustrates an alternate embodiment of a modular 
shuffling result and the fixed length data loss associated with the 
modular shuffling. 

Figure 15 illustrates one embodiment of variable length data 
buffering in a frame set. 

Figure 16 illustrates one embodiment of inter segment VL-data 
shuffling in accordance with the teachings of the present invention. 

Figure 17 is a flow diagram generally illustrating one 
embodiment of the data recovery process of the present invention. 

Figure 18 is a flow diagram of one embodiment of the Qbit and 
Motion Flag recovery process of the present invention. 

Figure 19 is a table illustrating one embodiment of candidate 
decodings. 

Figures 20a, 20b, 20c, 20d illustrate embodiments of 
measurements utilized in the Qbit and Motion Flag recovery process 
of Figure 18. 
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Figure 21 illustrates one embodiment of a table used to 
determine a square error probability function utilized in the Qbit and 
Motion Flag recovery process of Figure 18, 

Figure 22 illustrates one embodiment of a Qbit, Motion Flag and 
auxiliary information recovery process in accordance with the 
teachings of the present invention. 

Figure 23 illustrates the use of a post-amble in one embodiment 
of a bidirectional Qbit and Motion Flag recovery process. 

Figures 24a, 24b and 24c illustrate an alternate embodiment for 
evaluating candidate decodings. 

Figure 25 illustrates the use of smoothness measures in 
accordance with the teachings of one embodiment of the present 
invention. 

Figures 26a, 26b, 26c, 26d and 26e illustrate an alternate 
embodiment of a process for evaluating candidate decodings. 

Figure 27a illustrates an alternate process for evaluating 
candidate decodings and Figure 27b illustrates one embodiment for 
determining weighting values. 

DETAILED DESCRIPTION 
The present invention provides a method for coding and 
arranging a signal stream to provide for a robust error recovery. In the 
following description, for purposes of explanation, numerous details 
are set forth, in order to provide a thorough understanding of the 
present invention. However, it will be apparent to one skilled in the 
art that these specific details are not required in order to practice the 
present invention. In other instances, well known electrical 
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structures and circuits are shown in block diagram form in order not 
to obscure the present invention unnecessarily. 

The signal processing methods and structures are described 
from the perspective of one embodiment in which the signals are 
video signals. However, it is contemplated that the methods and 
apparatus described herein are applicable to a variety of types of signals 
including audio signals or other digital bitstreams of data, wherein 
each signal is composed of multiple signal elements. Furthermore the 
embodiment of the process described herein utilizes the Adaptive 
Dynamic Range Coding ("ADRC") process to compress data; however 
a variety of coding techniques and algorithms may be used. For a 
more detailed discussion on ADRC, see "Adaptive Dynamic Range 
Coding Scheme for Future HDTV Digital VTR", Kondo, Fujimori and 
Nakaya, Fourth International Workshop on HDTV and Beyond, 
September 4-6, 1991, Turin, Italy. 

In the above paper, three different kinds of ADRC are 
explained. These are achieved according to the following equations: 



Non-edge-matching ADRC: 

DR = MAX - MIN + 1 
(x-MIN + 0.5) 2 



<7 = 



x = 



Q 



DR 



(q + 0.5)«PR 



2 Q 



MIN 



Edge-matching ADRC: , 



DR = MAX - MIN 

\x-MlN)*{2 Q - 
DR 



x = 



S-2R + min + 

2 Q -1 



0.5 J 
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Multi-stage ADRC: 




DR = MAX - MIN + 1 

_\{x- MIN + 0.5) 2 Q 



Where MAX' is the averaged value of in the case of q = 2 Q - 1 ; 
MIN' is the averaged value of x in the case of q = 0; and 

DR* = MAX" - AfflV' 



where MAX represents the maximum level of a block, MIN 
represents the minimum level of a block, x represents the signal level 
of each sample, Q represents the number of quantization bits, q 
represents the quantization code (encoded data), x represents the 
decoded level of each sample, and the square brackets [ ] represent a 
truncation operation performed on the value within the square 
brackets. 

The signal encoding, transmission, and subsequent decoding 
processes are generally illustrated in Figure 1. Signal 100 is a data 
stream input to Encoder 110. Encoder 110 follows the Adaptive 
Dynamic Range Coding ("ADRC") compression algorithm and 
generates Packets 1, . . . N for transmission along Transmission Media 
135. Decoder 120 receives Packets 1, . . . N from Transmission Media 
135 and generates Signal 130. Signal 130 is a reconstruction of Signal 
100. 



(x-MIN') (2 Q -\) 
DR' 



+ 0.5 



x 



q PR' 



•f MIN' + 0.5 
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Encoder 110 and Decoder 120 can be implemented a variety of 
ways to perform the functionality described herein. In one 
embodiment, Encoder 110 and/or Decoder 120 are embodied as 
software stored on media and executed by a general purpose or 
specifically configured computer system, typically including a central 
processing unit, memory and one or more input/output devices and 
co-processors. Alternately, the Encoder 110 and /or Decoder 120 may be 
implemented as logic to perform the functionality described herein. 
In addition, Encoder 110 and/or Decoder 120 can be implemented as a 
combination of hardware, software or firmware. 

In the present embodiment Signal 100 is a color video image 
comprising a sequence of video frames, each frame including 
information representative of an image in an interlaced video system. 
Each frame is composed of two fields, wherein one field contains data 
of the even lines of the image and the other field containing the odd 
lines of the image. The data includes pixel values which describe the 
color components of a corresponding location in the image. For 
example, in the present embodiment, the color components consist of 
the luminance signal Y, and color difference signals U, and V. It is 
readily apparent the process of the present invention can be applied to 
signals other than interlaced video signals. Furthermore, it is 
apparent that the present invention is not limited to implementations 
in the Y, U, V color space, but can be applied to images represented in 
other color spaces. 

Referring back to Figure 1, Encoder 110 divides the Y, U, and V 
signals and processes each group of signals independently in 
accordance with the ADRC algorithm. The following description, for 
purposes of simplifying the discussion, describes the processing of the 
Y signal; however, the encoding steps are replicated for the U and V 
signals. 

In the present embodiment, Encoder 110 groups Y signals across 
two subsequent frames, referred to herein as a frame pair, of Signal 100 
into three dimensional blocks ("3D") blocks. For one embodiment, a. 
3D block is generated from grouping two 2D blocks from the same 
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localized area across a given frame pair, wherein a two dimensional 
2D block is created by grouping localized pixels within a frame or a 
field. It is contemplated that the process described herein can be 
applied to different block structures. The grouping of signals will be 
further described in the image-to-block mapping section below. 

Continuing with the present embodiment, for a given 3D block, 
Encoder 110 calculates whether there is a change in pixel values 
between the 2D blocks forming the 3D block. A Motion Flag is set if 
there are substantial changes in values. As is known in the art, use of 
a Motion Flag allows Encoder 110 to reduce the number of 
quantization codes when there is localized image repetition within 
each frame pair. Encoder 110 also detects the maximum pixel 
intensity value ("MAX") and the minimum pixel intensity value 
("MINT) within a 3D block. Using values MAX and MIN, Encoder 110 
calculates the dynamic range ("DR") for a given 3D block of data. For 
one embodiment DR ~ MAX - MIN + 1 in the case of non-edge- 
matching ADRC. For edge-matching ADRC, DR = MAX - MIN. In an 
alternative embodiment, Encoder 110 encodes signals on a frame by 
frame basis for a stream of frames representing a sequence of video 
frames. In another embodiment, Encoder 110 encodes signals on a 
field by field basis for a stream of fields representing a sequence of 
video fields. Accordingly, Motion Flags are not used and 2D blocks 
are used to calculate the MIN, MAX, and DR values. 

In the present embodiment, Encoder 110 references the 
calculated DR against a threshold table (not shown) to determine the 
number of quantization bits ("Qbits") used to encode pixels within the 
block corresponding to the DR. Encoding of a pixel results in a 
quantization code ( M Q code"). The Q codes are the relevant 
compressed image data used for storage or transmission purposes. 

In one embodiment, the Qbit selection is derived from the DR 
of a 3D block. Accordingly, all pixels within a given 3D block are 
encoded using the same Qbit, resulting in a 3D encoded block. The 
collection of Q codes, MIN, Motion Flag, and DR for a 3D encoded 
block is referred to as a 3D ADRC block. Alternately, 2D blocks are 
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encoded and the collection of Q codes, MEN, and DR for a given 2D 
block results in 2D ADRC blocks. 

A number of threshold tables can be implemented. In one 
embodiment, the threshold table consists of a row of DR threshold 
values. A Qbit corresponds to the number of quantization bits used to 
encode a range of DR values between two adjacent DRs within a row 
of the threshold table. In an alternative embodiment, the threshold 
table includes multiple rows and selection of a row depends on the 
desired transmission rate. Each row in the threshold table is identified 
by a threshold index. A detailed description of one embodiment of 
threshold selection is described below in the discussion of partial 
buffering. A further description of ADRC encoding and buffering is 
disclosed in US Patent no. 4,722,003 entitled "High Efficiency Coding 
Apparatus" and US Patent no. 4,845,560 also entitled "High Efficiency 
Coding Apparatus", assigned to the assignee of the present invention. 

Here forth the Q codes are referred to as variable length data 
("VL-data"). In addition, the DR, MIN, and Motion Flag are referred 
to as block attributes. The block attributes, together with the threshold 
index, constitute the fixed length data ("FL-data"). Furthermore, in 
view of the above discussion, the term block attribute describes a 
parameter associated with a component of a signal element, wherein a 
signal element includes multiple components. 

In an alternate embodiment, the FL-data includes a Qbit code. The 
advantage is that the Qbit information does not have to be derived 
from the DR during the decoding process. Thus, if the DR 
information is lost or damaged, the Qbit information can still be 
determined from the Qbit code. Furthermore, if the Qbit code is lost 
or damaged, the Qbit information can be derived from DR. Thus the 
requirement to recover the DR and Qbit is reduced. 

The disadvantage to including the Qbit code is the additional 
bits to be transmitted for each ADRC block. However, in one 
embodiment, Qbit codes for groups of ADRC blocks are combined, for 
example, in accordance with a function such as addition or 
concatenation. For example, if ADRC blocks are grouped in threes and 
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if the Qbit values for each ADRC block are respectively 3, 4 and 4, the 
summed value that is encoded into the FL-data is 11. Thus the 
number of bits required to represent the sum is less than the number 
of bits required to represent each individual value and undamaged 
Qbit values of the group can be used to determine the Qbit value 
without performing a Qbit recovery process such as the one described 
subsequently. 

Other embodiments are also contemplated. For example, 
Motion Flag data may also be encoded. A tag with Qbit and Motion 
Flag data can be generated and used to reference a table of codes. The 
configuration and function of the coding can vary according to 
application. 

Frames, block attributes, and VL-data describe a variety of 
components within a video signal. The boundaries, location, and 
quantity of these components are dependent on the transmission and 
compression properties of a video signal. In the present embodiment, 
these components are varied and shuffled within a bitstream of the 
video signal to ensure a robust error recovery during transmission 
losses. 

For illustrative purposes, the following description provides for 
a 1/6 consecutive packet transmission loss tolerance, pursuant to an 
ADRC encoding and shuffling of a video signal. Accordingly, the 
following definition and division of components exist for one 
embodiment. Other embodiments also are contemplated. A data set 
includes a partition of data of a video or other type of data signal. 
Thus, in one embodiment, a frame set is a type of data set that 
includes one or more consecutive frames. A segment includes a 
memory with the capacity to store a one-sixth division of the Q codes 
and block attributes included in a frame set. Further, a buffer includes 
a memory with the capacity to store a one-sixtieth division of the Q 
codes and block attributes included in a frame set. The shuffling of 
data is performed by interchanging components within segments 
and/ or buffers. Subsequently, the data stored in a segment is used to 
generate packets of data for transmission. Thus, in the following 
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description if a segment is lost all the packets generated from the 
segment are lost during transmission. Similarly, if a fraction of a 
segment is lost then a corresponding number of packets generated 
from the segment are lost during transmission. 

Although, the following description refers to a 1/6 consecutive 
packet loss for data encoded using ADRC encoding, it is contemplated 
that the methods and apparatus described herein are applicable to a 
design of a 1/n consecutive packets loss tolerance coupled to a variety 
of encoding/decoding schemes. 

Figure 2 illustrates one embodiment of Packet Structure 200 
used for the transmission of the data across point-to-point connections 
as well as networks. Packet Structure 200 is generated by Encoder 110 
and is transmitted across Transmission Media 135. For one 
embodiment, Packet Structure 200 comprises five bytes of header 
information, eight DR bits, eight MIN bits, a Motion Flag bit, a five bit 
threshold index, and 354 bits of Q codes. The packet structure 
described herein is illustrative and may typically be implemented for 
transmission in an asynchronous transfer mode ("ATM") network. 
However, the present invention is not limited to the packet structure 
described and a variety of packet structures that are used in a variety of 
networks can be utilized. 

As noted earlier, Transmission Media (e.g., media) 135 is not 
assumed to provide error-free transmission and therefore packets may 
be lost or damaged. As noted earlier, conventional methods exist for 
detecting such loss or damage, but substantial image degradation will 
generally occur. The system and methods of the present invention 
therefore teach source coding to provide robust recovery from such 
loss or damage. It is assumed throughout the following discussion 
that a burst loss, that is the loss of several consecutive packets, is the 
most probable form of error, but some random packet losses might 
also occur. 

To ensure a robust recovery for the loss of one or more 
consecutive packets of data, the system and methods of the present - 
invention provide multiple level shuffling. In particular, the FL-data 
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and the VL-data included in a transmitted packet comprise data from 
spatially and temporally disjointed locations of an image. Shuffling 
data ensures that any burst error is scattered and facilitates error 
recovery. As will be described below, the shuffling allows recovery of 
block attributes and Qbit values. 

Data Encoding/Decoding 
Figure 3 is a flow diagram illustrating one embodiment of the 
encoding process performed by Encoder 110. Figure 3 further describes 
an overview of the shuffling process used to ensure against image 
degradation and to facilitate a robust error recovery. 

In step one of Figure 3, an input frame set, also referred to as a 
display component, is decimated to reduce the transmission 
requirements. The Y signal is decimated horizontally to three- 
quarters of its original width and the U and V signals are each 
decimated to one-half of their original height and one-half of their 
original width. This results in a 3:1:0 video format with 3960 Y blocks, 
660 U blocks and 660 V blocks in each frame pair. As noted earlier, the 
discussion will describe the processing of Y signals; however, the ' 
process is applicable to the U and V signals. At step two, the two Y 
frame images are mapped to 3D blocks. At step three, 3D blocks are 
shuffled. At step four, ADRC buffering and encoding is used. At step 
five, encoded Y, U and V blocks are shuffled within a buffer. 

At step six, the VL-data for a group of encoded 3D blocks and 
their corresponding block attributes are shuffled. At step seven, the 
FL-data is shuffled across different segments. At step eight, post-amble 
filling is performed in which variable space at the end of a buffer is 
filled with a predetermined bitstream. At step nine, the VL-data is 
shuffled across different segments. 

For illustrative purposes the following shuffling description 
provides a method for manipulation of pixel data before and after 
encoding. For an alternative embodiment, independent data values 
are shuffled /deshuf fled via hardware. In particular, the hardware 
maps the address of block values to different addresses to implement 
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the shuffling/ deshuf fling process. However, address mapping is not 
possible for data dependent values because shuffling has to follow the 
processing of data. The intra group VL-data shuffling described below 
includes the data dependent values. Further, for illustrative purposes 
the following shuffling description occurs on discrete sets of data. 
However, for alternative embodiments a signal is defined based on 
multiple data levels ranging from bits, to pixels, and to frames. 
Shuffling is possible for each level defined in the signal and across 
different data levels of the signal. 

Figure 4 is a flow diagram illustrating one embodiment of 
decoding process performed by Decoder 120. Preferably, the 
conversion and de-shuffling processes are the inverse of the processes 
represented in Figure 3. Figure 4, further describes, in different 
combinations of Qbit, Motion Flag, DR, MIN and pixel data, an 
innovative process for error recovery. The error recovery process is 
described below in different combinations for different embodiments, 
Qbit, Motion Flag, DR, MIN and pixel recovery, 

Image-to-Block Mapping 

In the present embodiment, a single frame typically comprises 
5280 2D blocks wherein each 2D block comprises 64 pixels. Thus, a 
frame pair comprises 5280 3D blocks as a 2D block from a first frame 
and a 2D block from a subsequent frame are collected to form a 3D 
block. 

Image-to-block mapping is performed for the purpose of 
dividing a frame or frame set of data into 2D blocks or 3D blocks 
respectively. Moreover, image-to-block mapping includes using a 
complementary and/or interlocking pattern to divide pixels in a 
frame to facilitate robust error recovery during transmission losses. 
However, to improve the probability that a given DR value is not too 
large, each 2D block is constructed from pixels in a localized area. 

Figure 5 illustrates one embodiment of an image-to-block 
mapping process for an exemplary 16 pixel section of an image. Image 
500 comprises 16 pixels forming a localized area of a single frame. 
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Each pixel in Image 500 is represented by an intensity value. For 
example, the pixel in the top left hand side of the image has an 
intensity value equal to 100 whereas the pixel in the bottom right 
hand side of the image has an intensity value of 10. 

In one embodiment, pixels from different areas of Image 500 are 
used to create 2D Blocks 510, 520, 530, and 540. 2D Blocks 510, 520, 530, 
and 540 are encoded, shuffled (as illustrated below), and transmitted. ' 
Subsequent to transmission, 2D Blocks 510, 520, 530, and 540 are 
recombined and used to form Image 550. Image 550 is a reconstruction 
of Image 500. 

To ensure accurate representation of Image 500 despite a 
possible transmission loss, Figure 5 is an interlocking complementary 
block structure, one embodiment of which is illustrated in Figure 5, is 
used to reconstruct Image 550. In particular, the pixel selection used 
to create 2D Blocks 510, 520, 530, and 540 ensures that a complementary 
and/or interlocking pattern is used to recombine the blocks when 
Image 550 is reconstructed. Accordingly, when a particular 2D block's 
attribute is lost during transmission, contiguous sections of Image 550 
are not distorted during reconstruction. For example, as illustrated in 
Figure 5 the DR of 2D Block 540 is lost during data transmission. 
However, during reconstruction of Image 550, the decoder utilizes 
multiple neighboring pixels of neighboring blocks through which a 
DR can be recovered for the missing DR of 2D Block 540. In addition, 
as will be subsequently described, the combination of complementary 
patterns and shifting increases the number of neighboring pixels, 
preferably maximizing the number of neighboring pixels that 
originate from other blocks, significantly improving DR and MIN 
recovery. 

Figure 5a illustrates one embodiment of a shuffling pattern' 
used to form 2D blocks in one embodiment of the image-to-block 
mapping process. An image is decomposed into two sub-images, Sub- 
Image 560 and Sub-Image 570, based on alternating pixels. Rectangular 
shapes are formed in Sub-Image 560 to delineate the 2D block 
boundaries. For purposes of discussion, the 2D blocks are numbered 
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0, 2, 4, 7, 9, 11, 12, 14, 16, 19, 21, and 23. Tile 565 illustrates the pixel 
distribution for a 2D block within Sub-Image 560. 

In Sub-Image 570, the 2D block assignment is shifted by eight 
pixels horizontally and four pixels vertically. This results in a wrap 
around 2D block assignment and overlap when Sub-Images 560 and 
570 are combined during reconstruction. The 2D blocks are numbered 

1, 3, 5, 6, 8, 10, 13, 15, 17, 18, 20, and 22. Tile 575 illustrates the pixel 
distribution for a 2D block within Sub-Image 570. Tile 575 is the 
complementary structure of Tile 565. Accordingly, when a particular 
block's attribute is lost during transmission, neighboring pixels 
through which a block attribute can be recovered for the missing 2D 
block exists. Additionally, an overlapping 2D block of pixels with a 
similar set of block attributes exist. Therefore, during reconstruction 
of the image the decoder has multiple neighboring pixels from 
adjacent 2D blocks through which a lost block attribute can be 
recovered. 

Figure 6 illustrates other complementary and interlocking 2D 
block structures. Other structures may also be utilized. Similar to 
Figure 5, these 2D block structures illustrated in Figure 6, ensure 
surrounding 2D blocks are present despite transmission losses for a 
given 2D block. However, Patterns 610a, 610b, and 610d use horizontal 
and /or vertical shifting during the mapping of pixels to subsequent 2D 
blocks. Horizontal shifting describes shifting the tile structure in the 
horizontal direction a predetermined number of pixels prior to 
beginning a new 2D block boundary. Vertical shifting describes 
shifting the tile structure in the vertical direction a predetermined 
number of pixels prior to beginning a new 2D block boundary. In 
application, horizontal shifting only may be applied, vertical shifting 
may only be applied, or a combination of horizontal and vertical 
shifting may be applied. 

Pattern 610a illustrates a spiral pattern used for image-to-block 
mapping. The spiral pattern follows a horizontal shifting to create 
subsequent 2D blocks during the image-to-block mapping process. 
Patterns 610b and 610d illustrate complementary patterns wherein 
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pixel selection is moved by a horizontal and vertical shifting to create 
subsequent 2D blocks during the image-to-block mapping process. 
Further, Patterns 610b and 610d illustrate alternating offsets on pixels 
selection between 2D blocks. Pattern 610c illustrates using an irregular 
sampling of pixels to create a 2D block for image-to-block mapping. 
Accordingly, the image-to-block mapping follows any mapping 
structure provided a pixel is mapped to a 2D block only once. 

Figure 5, Figure 5a and Figure 6 describe image-to-block 
mapping for 2D block generation. It is readily apparent that the 
processes are applicable to 3D blocks. As described above, 3D block 
generation follows the same boundary definition as a 2D block, 
however the boundary division extends across a subsequent frame 
resulting in a 3D block. In particular, a 3D block is created by collecting 
the pixels used to define a 2D block in a first frame together with pixels 
from a 2D block in a subsequent frame. In one embodiment, both 
pixels in the 2D block from the first frame and the 2D block from the 
subsequent frame are from the exact same location. 

Intra Frame Set Block Shuffling 

The pixels values for a given image are closely related for a 
localized area. However, in another area of the same images the pixel 
values may have significantly different values. Thus, subsequent to 
encoding the DR and MIN values for spatially close 2D or 3D blocks in 
a section of an image have similar values, whereas the DR and MIN 
values for blocks in another section of the image may be significantly 
different. Accordingly, when buffers are sequentially filled with 
encoded data from spatially close 2D or 3D blocks of an image, a 
disproportionate usage of buffer space occurs. Intra frame set block 
shuffling occurs prior to ADRC encoding and includes shuffling the 
2D or 3D blocks generated during the image-to-block mapping process. 
This shuffling process ensures an equalized buffer usage during a 
subsequent ADRC encoding. 

Figures 7a - 7d illustrate one embodiment of shuffling 3D Y- 
blocks. The 3D Y-blocks in Figures 7a-7d are generated from applying 
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the image-to-block mapping process described above to a frame pair 
containing only Y signals. The 3D Y-blocks are shuffled to ensure that 
the buffers used to store the encoded frame pair contain 3D Y-blocks 
from different parts of the frame pair. This leads to similar DR 
distribution during ADRC encoding. A similar DR distribution 
within each buffer leads to consistent buffer utilization. 

Figure 7a -7d also illustrate 3D block shuffling using physically 
disjointed 3D blocks to ensure that transmission loss of consecutive 
packets results in damaged block attributes scattered across the image, 
as opposed to a localized area of the image. 

The block shuffling is designed to widely distribute block 
attributes in the event of small, medium, or large, burst packet losses 
occur. In the present embodiment, a small burst loss is thought of as 
one where a few packets are lost; a medium loss is one in which the 
amount of data that can be held in one buffer is lost; and a large loss is 
one in which the amount of data that can be held in one segment is 
lost. During the 3D block shuffling each group of three adjacent blocks 
are selected from relatively remote parts of the image. Accordingly, 
during the subsequent intra group VL-data shuffling (to be detailed 
later), each group is formed from 3D blocks that have differing 
statistical characteristics. Distributed block attribute losses allow for a 
robust error recovery because a damaged 3D block is surrounded by 
undamaged 3D blocks and the undamaged 3D blocks can be used to 
recover lost data. 

Figure 7a illustrates a frame pair containing 66 3D Y-blocks in 
the horizontal direction and 60 3D Y-blocks in the vertical direction. 
The 3D Y-blocks are allocated into Segments 0-5. As illustrated, the 
3D Y-block assignment follows a two by three column section such 
that one 3D Y-block from each section is associated with a segment. 
Thus, if no further shuffling is performed and a burst loss of the first 
880 packets occurs, all the block attributes associated with Segment 0 
are lost. However, as later described, FL-data shuffling is performed to 
further disperse block attribute losses. 
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Figure 7b illustrates the scanning order of 3D Y-blocks 
numbered "0" used to enter into Segment 0. Each "0" 3D Y-block of 
Figure 7a is numbered 0, 1, 2, 3, . . . ., 659 to illustrate their location in 
the stream that is inputted into Segment 0. Using the block 
numbering to allocate segment assignment the remaining 3D Y-blocks 
are inputted into Segments 1-5, thus resulting in a frame pair 
shuffled across multiple segments. 

Figure 7c illustrates the 660 3D Y-blocks comprising one 
segment. The 3D Y-blocks numbered 0-65 are inputted into Buffer 0. 
Similarly the 3D Y-blocks adjacent to the numbered 3D Y-blocks are 
inputted into Buffer 1. The process is repeated to fill Buffers 2-9. 
Accordingly, damage to a buffer during data transmission results in 
missing 3D Y-blocks from different parts of the image. 

Figure 7d illustrates the final ordering of the "0" 3D Y-blocks 
across a buffer. 3D Y-blocks 0, 1, and 2 occupy the first three positions 
in the buffer. The process is repeated for the rest of the buffer. 
Accordingly, loss of three 3D Y-blocks during data transmission results 
in missing 3D Y-blocks from distant locations within the image. 

Figures 7a-d illustrate one embodiment of 3D block 
distributions for 3D Y-blocks of a frame set. In alternative 
embodiments, however, 3D block distributions for 3D U-blocks and 3D 
V-blocks are available. The 3D U-blocks are generated from applying 
the image-to-block mapping process, described above, to a frame set 
containing only U signals. Similarly, 3D V-blocks are generated from 
applying the image-to-block mapping process to a frame set containing 
only V signals. Both the 3D U-block and the 3D V-block follow the 3D 
Y-block distribution described above. However, as previously 
described, the number of 3D U-blocks and 3D V-blocks each have a 1:6 
proportion to 3D Y-blocks. 

Figures 7a-d are used to illustrate one embodiment of intra 
frame set block shuffling for a Y signal such that burst error of up to 
1/6 of the packets lost during transmission is tolerated and further 
ensures equalized buffer use. It will be appreciated by one skilled in . 
the art that segment, buffer, and ADRC block assignments can be 
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varied to ensure against 1 /n burst error loss or to modify buffer 
utilization. 

Partial Buffering 

As illustrated in Figure 3, the ADRC encoding and buffering 
processes occur in step four. Dependent on the encoding technique, 
2D or 3D blocks generated during the image-to-block mapping process 
are encoded resulting in 2D or 3D ADRC blocks. A 3D ADRC block, 
contains Q codes, a MIN value, a Motion Flag, and a DR. Similarly, a 
2D ADRC block contains Q codes, a MIN, and a DR. A 2D ADRC block, 
however, does not include a Motion Flag because the encoding is 
performed on a single frame or a single field. 

A number of buffering techniques are found in the prior art (see 
for example, High Efficiency Coding Apparatus, U.S. Patent 4,845,560 
of Kondo et. al. and High Efficiency Coding Apparatus, U.S. Patent 
4,722,003 of Kondo). Both High Efficiency Coding Apparatus patents 
are hereby incorporated by reference. 

The partial buffering process set forth below, describes an 
innovative method for determining the encoding bits used in ADRC 
encoding. In particular, partial buffering describes a method of 
selecting threshold values from a threshold table designed to provide 
a constant transmission rate between remote terminals while 
restricting error propagation. In an alternative embodiment, the 
threshold table is further designed to provide maximum buffer 
utilization. In one embodiment, a buffer is a memory that stores a 
one-sixtieth division of encoded data from a given frame set. The 
threshold values are used to determine the number of Qbits used to 
encode the pixels in 2D or 3D blocks generated from the image-to- 
block mapping process previously described. 

The threshold table includes rows of threshold values, also 
referred to as a threshold set, and each row in the threshold table is 
indexed by a threshold index. In one embodiment, the threshold table 
is organized with threshold sets that generate a higher number of Q . 
code bits located in the upper rows of the threshold table. Accordingly, 
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for a given buffer having a predetermined number of bits available, 
Encoder 110 moves down the threshold table until a threshold set that 
generates less than a predetermined number of bits is encountered. 
The appropriate threshold values are used to encode the pixel data in 
the buffer. 

In one embodiment, a transmission rate of no more than 30 
Mbps is desired. The desired transmission rate results in 31,152 bits 
available for VL-data storage in any given buffer. Accordingly, for 
each buffer a cumulative DR distribution is computed and a threshold 
set is selected from the threshold table to encode the pixels in 3D or 2D 
blocks into VL-data. 

Figure 8 illustrates one embodiment of selected threshold 
values and the DR distribution for Buffer 0. The vertical axis of Figure 
8 includes the cumulative DR distribution. For example, the value 
"b" is equal to the number of 3D or 2D blocks whose DR is greater 
than or equal to L3. The horizontal axis includes the possible DR 
values. In one embodiment, DR values range from 0 to 255. 
Threshold values L4, L3, L2, and Li describe a threshold set used to 
determine the encoding of a buffer. 

In one embodiment, all blocks stored in Buffer 0 are encoded 
using threshold values L4, L3, L2, and Li- Accordingly, blocks with DR 
values greater than L4 have their pixel values encoded using four bits. 
Similarly, all pixels belonging to blocks with DR values between L3 
and L4 are encoded using three bits. All pixels belonging to blocks 
with DR values between L2 and L3 are encoded using two bits. All 
pixels belonging to blocks with DR values between I4 and L2 are 
encoded using one bit. Finally, all pixels belonging to blocks with DR 
values smaller than Li are encoded using zero bits. L4, L3, L2/ and Li 
are selected such that the total number of bits used to encode all the 
blocks in Buffer 0 is as close as possible to a limit of 31,152 bits without 
exceeding the limit of 31,152. 

Figure 8a illustrates the use of partial buffering in one 
embodiment. Frame 800 is encoded and stored in Buffers 0-59. 
Provided a transmission error inhibits data recovery, the decoding 
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process is stalled for Frame 800 until error recovery is performed on 
the lost data. However, partial buffering restricts the error 
propagation~within a buffer, thus allowing decoding of the remaining 
buffers. In one embodiment, a transmission error inhibits the Qbit 
and Motion Flag recovery for Block 80 in Buffer 0. Partial buffering 
limits the error propagation to the remaining blocks within Buffer 0. 
Error propagation is limited to Buffer 0 because the end of Buffer 0 and 
the beginning of Buffer 1 are known due to the fixed buffer length. 
Accordingly, Decoder 120 can begin processing of blocks within Buffer 
1 without delay. Additionally, the use of different threshold sets to 
encode different buffers allows Encoder 110 to maximize/control the 
number of Q codes bits included in a given buffer, thus allowing a 
higher compression ratio. Furthermore, the partial buffering process 
allows for a constant transmission rate because Buffers 0-59 consist of 
a fixed length. 

In one embodiment, a buffer's variable space is not completely 
filled with Q code bits because a limited number of threshold sets exist. 
Accordingly, the remaining bits in the fixed length buffer are filled 
with a predetermined bitstream pattern referred to as a post-amble. As 
will be described subsequently, the post-amble enables bidirectional 
data recovery because the post-amble delineates the end of the VL-data 
prior to the end of the buffer. 

Intra Buffer YUV Block Shuffling 

Y, U, and V, signals each have unique statistical properties. To 
improve the Qbit and Motion Flag recovery process (described below) 
the Y, U, and V signals are multiplexed within a buffer. Accordingly, 
transmission loss does not have a substantial effect on a specific signal. 

Figure 9 illustrates one embodiment of the intra buffer YUV 
block shuffling process in which YUV ADRC blocks are derived from 
the Y, U, and V signals respectively. Buffer 900 illustrates the ADRC 
block assignments after intra frame set block shuffling. Buffer 900 
comprises 66 Y-ADRC blocks followed by 11 U-ADRC blocks which are 
in turn followed by 11 V-ADRC blocks. Buffer 910 shows the YUV 
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ADRC block organization after intra buffer YUV block shuffling. As 
illustrated, three Y-ADRC blocks are followed by a U-ADRC block or 
three Y-ADRC blocks are followed by a V-ADRC block. Intra buffer 
YUV block shuffling reduces similarity between adjacent block's 
bitstreams within the buffer. Alternative embodiments of intra buffer 
YUV block shuffling with a different signal, i.e., YUV ratios or other 
color spaces are possible dependent on the initial image format. 

Intra Group VL-Data Shuffling 

Intra group VL-data shuffling comprises three processing steps. 
The three processing steps include Q code concatenation, Q code 
reassignment, and randomizing concatenated Q codes. Figure 10 
illustrates one embodiment of intra group VL-data shuffling wherein 
three processing steps are applied consecutively to Q codes stored in a 
buffer. In alternative embodiments, a subset of the processing step is 
applied in intra group VL-data shuffling. Each processing step 
independently assists in the error recovery of data lost during 
transmission. Accordingly, each processing step is described 
independently. A detailed description of error recovery is provided 
below in the discussion of data recovery. 

1. Q code concatenation 

Q code concatenation ensures that groups of ADRC blocks are 
decoded together. Group decoding facilitates error recovery because 
additional information is available from neighboring blocks during 
the data recovery process detailed below. For one embodiment, Q code 
concatenation is applied independently to each group of three ADRC 
blocks stored in a buffer. In an alternative embodiment, a group 
includes ADRC block(s) from different buffers. The concatenation of 
Q codes across three ADRC blocks is described as generating one 
concatenated ADRC tile. Figure 11 and Figure 11a illustrate one 
embodiment of generating concatenated ADRC tiles. 

Figure 11 illustrates one embodiment of generating a 
concatenated ADRC tile from 2D ADRC blocks. Specifically, the 
concatenation is performed for each Q code (qo - q63) included in 2D 
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ADRC Blocks 0, 1, and 2 resulting in the sixty four Q codes of 
Concatenated ADRC Tile A. For example, the first Q code qo,o (0th 
quantized value) of 2D ADRC Block 0 is concatenated to the first Q 
code qo,i of 2D ADRC Block 1. The two concatenated Q codes are in 
turn concatenated to the first Q code qo,2 of 2D ADRC Block 2, thus 
resulting in Qo of Concatenated ADRC Tile A. The processes is 
repeated until Qe^ is generated. Alternatively, the generation of Qi in 
Concatenated ADRC Tile A is described by the equation 

Qi = [ qi,0 / qu / qi,2 ] i= 0, 1, 2 63 

Additionally, associated with each Qi in Concatenated ADRC 
Tile A there is a corresponding number of N bits that represents the 
total number of bits concatenated to generate a single Qi. 

Figure 11a illustrates one embodiment of generating a 
concatenated ADRC tile from frame pairs including motion blocks. A 
motion block is a 3D ADRC block with a set Motion Flag. The Motion 
Flag is set when a predetermined number of pixels within two 2D 
blocks structure created by image-to-block mapping process described 
earlier, change in value between a first frame and a subsequent frame. 
In an alternative embodiment, the Motion Flag is set when the 
maximum value of each pixel change between the 2D block of a first 
frame and a subsequent frame exceeds a predetermined value. In 
contrast, non-motion (i.e., stationary) block includes a 3D ADRC block 
with a Motion Flag that is not set. The Motion Flag remains un-set 
when a predetermined number of pixels within the two 2D blocks of a 
first frame and a subsequent frame do not change in value. In an 
alternative embodiment, the Motion Flag remains un-set when the 
maximum value of each pixel change between a first frame and a 
subsequent frame does not exceed a predetermined value. 

A motion block includes Q codes from an encoded 2D block in a 
first frame and an encoded 2D block in a subsequent frame. The 
collection of Q codes corresponding to a single encoded 2D block are 
referred to as an ADRC tile. Accordingly, a motion block generates 
two ADRC tiles. However, due to the lack of motion, a stationary 
block need only include one-half of the number of Q codes of a 
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motion block, thus generating only one ADRC tile. In the present 
embodiment, the Q codes of a stationary block are generated by 
averaging corresponding pixels values between a 2D block in a first 
frame and a corresponding 2D block in a subsequent frame. Each 
averaged pixel value is subsequently encoded resulting in the 
collection of Q codes forming a single ADRC tile. Accordingly, Motion 
Blocks 1110 and 1130 generate ADRC Tiles 0, 1, 3, and 4. Stationary 
Block 1120 generates ADRC Tile 2. 

The concatenated ADRC tile generation of Figure 11a 
concatenates the Q codes for ADRC Tiles 0-4 into Concatenated 
ADRC Tile B. Specifically, the concatenation is performed for each Q 
code (qo - q63> included in ADRC Tiles 0, 1, 2, 3 and 4 resulting in the 
sixty four Q codes of Concatenated ADRC Tile B. Alternatively, the 
generation of each Q code, Qi, in Concatenated ADRC Tile B is 
described by the mathematical equation 

Qi = [ qi,o / qu / qu , qu / qi,4 ] i= o, l, 2, . . . 63 

2. Q code reassignment 

Q code reassignment ensures that bit errors caused by 
transmission losses are localized within spatially disjointed pixels. In 
particular, during Q code reassignment, Q codes are redistributed and 
the bits of the redistributed Q codes are shuffled. Accordingly, Q code 
reassignment facilitates error recovery because undamaged pixels 
surround each damaged pixel. Furthermore, DR and MESJ recovery is 
aided because pixel damage is distributed evenly throughout an ADRC 
block, DR and MIN recovery is detailed below in the data recovery 
discussion. 

Figure 12 illustrates one embodiment of pixel corruption 
during the transmission loss of a 1/6 burst error loss. In particular, 2D 
ADRC Blocks 1210, 1220, and 1230 each include sixty four pixels 
encoded using three bits. Accordingly, each pixel, Po through P63/ of a 
2D ADRC block is represented by three bits. 2D ADRC Block 1210 
shows the bit loss pattern, indicated by a darkened square, of bits when 
the first bit of every six bits are lost. Similarly, the bit loss pattern 
when the second bit or fourth bit of every six bits are lost are shown in 
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2D ADRC Blocks 1220 and 1230, respectively. Figure 12 illustrates that 
without Q code reassignment one-half of all the pixels 2D ADRC 
Blocks 1210, 1220, and 1230 are corrupted for a 1/6 burst error loss. 

For one embodiment, Q code reassignment is applied 
independently to each concatenated ADRC tile stored in a buffer, thus 
ensuring that bit errors are localized within spatially disjointed pixels 
upon deshuffling. In an alternative embodiment, Q code 
reassignment is applied to each ADRC block stored in a buffer. 

Figure 12a illustrates one embodiment of Q code reassignment 
that generates a bitstream of shuffled Q code bits from a concatenated 
ADRC tile. Table 122 and Table 132 illustrate the Q code 
redistribution. Bitstreams 130 and 140 illustrate the shuffling of Q 
code bits. 

Table 122 shows the concatenated Q codes for Concatenated 
ADRC Tile A. Qo is the first concatenated Q code and Qe3 is the final 
concatenated Q code. Table 132 illustrates the redistribution of Q 
codes. For one embodiment Qo, Q6/ Ql2, Qi8, Q24, Q30/ Q36/ Q42/ Q48/ 
Q54, and Qeo are included in a first set, partition 0. Following Table 
132, the following eleven concatenated Q codes are included in 
partition 1. The steps are repeated for partitions 2 - 5. The boundary of 
a partition is delineated by a vertical line in Table 132. This disjointed 
spatial assignment of concatenated Q codes to six partitions ensures 
that a 1/6 burst error loss results in a bit loss pattern distributed across 
a group of consecutive pixels. 

Figure 12b illustrates one embodiment of the bit pattern loss 
created by the 1/6 burst error loss of redistributed Q codes. In 
particular, 2D ADRC blocks 1215, 1225, and 1235 each include sixty 
four pixels encoded using three bits. Accordingly, each pixel Po 
through P63, of each 2D ADRC block, is represented by three bits. In 2D 
ADRC Blocks 1215, 1225, and 1235 the bit loss pattern, indicated by a 
darkened square, is localized across a group of consecutive pixels. 
Accordingly, only eleven consecutive pixels within each 2D ADRC 
Block 1215, 1225, and 1235 are corrupted for a given segment loss. In 
an alternative embodiment, Q code assignment to partitions include Q 
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codes from different motion blocks, thus providing both a disjointed 
spatial and temporal assignment of Q codes to six segments. This 
results in additional undamaged spatial-temporal pixels during a 1/6 
burst error loss and further facilitates a more robust error recovery. 

Referring to Figure 12a, the bits of the redistributed Q codes in 
Table 132 are shuffled across a generated bitstream so that adjacent bits 
in the bitstream are from adjacent partitions. The Q code bits for all 
the partitions in Table 132 are concatenated into Bitstream 130. For a 
given partition adjacent bits in Bitstream 130 are scattered to every 
sixth bit location in the generated Bitstream 140. Accordingly, bits 
number zero through five, of Bitstream 140, include the first bit from 
the first Q code in each partition. Similarly, bits number six through 
eleven, of Bitstream 140, include the second bit from the first Q code 
in each partition. The process is repeated for all Q code bits. 
Accordingly, a 1/6 burst error loss will result in a spatially disjointed 
pixel loss. 

Figure 12c illustrates one embodiment of the bit pattern loss 
created by the 1/6 burst error loss of reassigned ( i.e. redistributed and 
shuffled) Q codes. In particular, 2D ADRC Blocks 1217, 1227, and 1237 
each include sixty four pixels encoded using three bits. Accordingly, 
each pixel Po through P63, of each 2D ADRC Block, is represented by 
three bits. In 2D ADRC Blocks 1217, 1227, and 1237, the bit loss pattern, 
indicated by a darkened square, is distributed across spatially disjointed 
pixels, thus facilitating pixel error recovery. 

3. Randomization of Q codes bits 

The Q code bits are randomized using a masking key to assist 
the decoder in recovering lost and damaged data. In particular, during 
encoding a key, denoted by KEY, is used to mask a bitstream of Q codes. 
Accordingly, the decoder must discern the correct values of KEY to 
unmask the bitstream of Q codes. 

In one embodiment, KEY is used to mask a bitstream of Q codes 
generated by the Q code reassignment of three ADRC blocks. As 
previously described, an ADRC block includes FL-data and Q codes. - 
Each key element ("di") of the masking key is generated by the 
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combination of the FL-data values and the number of quantization 
bits ("qi") associated with a corresponding ADRC block. In one 
embodiment, Motion Flags and Qbits are used to define a key. 
Accordingly, in this embodiment, the value of a key element is 
generated from the mathematical equation 

dj = 5- mi + qi where i = 0, 1, 2 and qi = 0, 1, 2, 3, 4 
The variable mj equals the Motion Flag. Accordingly, when the 
corresponding ADRC block is a stationary block, mj equals 0 and when 
the corresponding ADRC block is a motion block, mi equals 1. 
Furthermore, the variable qi represents the quantization bits used to 
encode the corresponding ADRC block. Accordingly, qi has a value of 
0, 1, 2, 3, or 4 for a four bit ADRC encoding technique. In one 
embodiment, KEY for a group of three ADRC blocks is defined with 
three key elements ("di") according to the following equation: 

KEY = d 0 + 10 • di + 100 • d 2 
Thus, during the recovery of Motion Flag or Qbit data possible key 
values are regenerated depending on the values used to create the 
masking keys. The regenerated key values are used to unmask the 
received bitstream of Q codes resulting in candidate decodings. A 
detailed description of regenerating key values and the selection of a 
specific candidate decoding is provided below in the discussion of data 
recovery. 

In an alternative embodiments, the masking key is generated 
form a variety of elements. Thus, providing the decoder with the 
specific information relating to an element without having to 
transmit the element across a transmission media. In one 
embodiment, DR or MIN values corresponding to an ADRC block are 
used to generate a masking key to mask the bitstream representing the 
ADRC block. 

Figures 10 - 12 illustrate intra group VL-data shuffling tolerated 
up to 1/6 packet data loss during transmission. It will be appreciated 
by one skilled in the art, that the number of total partitions and bit 
separation can be varied to ensure against 1/n burst error loss. 
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Inter Segment FL-Data Shuffling 

Inter segment FL-data shuffling describes rearranging block 
attributes among different segments. Rearranging block attributes 
provides for a distributed loss of data. In particular, when FL-data 
from a segment is lost during transmission the DR value, MIN value, 
and Motion Flag value lost do not belong to the same block. Figures 
13 and 14 illustrate one embodiment of inter segment FL-data 
shuffling. 

Figure 13 illustrates the contents of Segments 0 to 5. For one 
embodiment, each segment comprises 880 DRs, 880 MINs, 880 Motion 
Flags, and VL-data corresponding to 660 Y-blocks, 110 U-blocks, and 
110 V-blocks. As illustrated in graph MIN Shuffling 1300, the MIN 
values for Segment 0 are moved to Segment 2, the MIN values for 
Segment 2 are moved to Segment 4, and the MIN values for Segment 
4 are moved to Segment 0. Additionally, the MIN values for 
Segment 1 are moved to Segment 3, the MIN values for Segment 3 are 
moved to Segment 5, and the Motion Flag values for Segment 5 are 
moved to Segment 1. 

Figure 13a illustrates Motion Flag shuffling. As illustrated, in 
graph Motion Flag Shuffling 1305, the Motion Flag values for 
Segment 0 are moved to Segment 4, the Motion Flag values for 
Segment 2 are moved to Segment 0, and the Motion Flag values for 
Segment 4 are moved to Segment 2. Additionally, the Motion Flag 
values for Segment 1 are moved to Segment 5, the Motion Flag values 
for Segment 3 are moved to Segment 1, and the Motion Flag values 
for Segment 5 are moved to Segment 3. Loss pattern 1310 illustrates 
the FL-data loss after Segment 0 is lost during transmission. 

For a specific block attribute, both Figure 13 and Figure 13a 
illustrate shuffling all instances of the specific block attribute between 
segments. For example, in Figure 13 the 880 MIN values from 
Segment 0 are collectively exchanged with the 880 MIN values in 
Segment 2. Similarly, in Figure 13a the 880 Motion Flags for Segment 
0 are collectively exchanged with the 880 Motion Flags in Segment 4. 
During a transmission loss of consecutive packets, this collective 
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shuffling of block attributes results in a disproportionate loss of a 
specific block attributes for a block group. In one embodiment, a block 
group includes three ADRC blocks. 

Figure 14 illustrates one embodiment of a modular three 
shuffling process for DR, MIN, and Motion Flag values. A modular 
three shuffling describes a shuffling pattern shared across three blocks 
(i.e., a block group) in three different segments. The shuffling pattern 
is repeated for all block groups within the three different segments. 
However, a different shuffling pattern is used for different block 
attributes. Accordingly, the modular three shuffling process 
distributes block attributes over all three segments. In particular, for a 
given block group a modular three shuffling ensures that only one 
instance of a specific block attribute is lost during the transmission loss 
of a segment. Thus, during the data recovery process, described below, 
a reduced number of candidate decodings are generated to recover data 
loss within a block. 

As illustrated in DR Modular Shuffle 1410, a segment stores 880 
DR values. Accordingly, the DR values are numbered 0 - 879 
dependent on the block from which a given DR value is derived. In a 
modular three shuffling the FL-data contents of three segments are 
shuffled. A count of 0 - 2 is used to identify each DR value in the 
three segments identified for a modular shuffling. Accordingly, DR's 
belonging to blocks numbered 0, 3, 6, 9 . . . belong to Count 0. 
Similarly, DR's belonging to blocks numbered 1, 4, 7, 10, . . . belong to 
Count 1 and DR's belonging to blocks numbered 2, 5, 8, 11 . . . belong to 
Count 2. Thus, for a given count the DR values associated with that 
count are shuffled across Segment 0, 2, and 4. Similarly, the DR 
values associated with the same count are shuffled across Segments 1, 
3, and 5. 

In DR Modular Shuffle 1410, the DR values belonging to Count 
0 are left un-shuffled. The DR values belonging to Count 1 are 
shuffled. In particular, the Count 1 DR values in Segment A are 
moved to Segment B, the Count 1 DR values in Segment B are moved 
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to Segment C, and the Count 1 DR values in Segment C are moved to 
Segment A. 

The DR values belonging to Count 2 are also shuffled. In 
particular, the Count 2 DR values in Segment A are moved to 
Segment C, the Count 2 DR values in Segment B are moved to 
Segment A, and the Count 2 DR values in Segment C are moved to 
Segment B. 

MIN Modular Shuffle 1420 illustrates one embodiment of a 
modular three block attribute shuffling process for MIN values. A 
segment includes 880 MIN values. In MIN Modular Shuffle 1420, the 
shuffling pattern used for Count 1 and Count 2 in DR Modular 
Shuffle 1410 are shifted to Count 0 and Count 1. In particular, the 
shuffling pattern used for Count 1 in DR Modular Shuffle 1410 is 
applied to Count 0. The shuffling pattern used for Count 2 in DR 
Modular Shuffle 1410 is applied to Count 1 and the MIN values 
belonging to Count 2 are left un-shuffled. 

Motion Flag Modular Shuffle 1430 illustrates one embodiment 
of a modular three block attribute shuffling process for Motion Flag 
values. A segment includes 880 Motion Flag values. In Motion Flag 
Modular Shuffle 1430, the shuffling pattern used for Count 1 and 
Count 2 in DR Modular Shuffle 1410 are shifted to Count 2 and Count 
0 respectively. In particular, the shuffling pattern used for Count 2 in 
DR Modular Shuffle 1410 is applied to Count 0. The shuffling pattern 
used for Count 1 in DR Modular Shuffle 1410 is applied to Count 2 
and the Motion Flag values belonging to Count 1 are left un-shuffled. 

Figure 14a illustrates the modular shuffling result of Modular 
Shuffles 1410, 1420, and 1430. Modular Shuffle Result 1416 shows each 
attribute destination of blocks belonging to Segment 0. In this 
example, Segment 0 corresponds to Segment A of Figure 14. This 
destination is defined according to Modular Shuffles 1410, 1420, and 
1430 of Figure 14. Figure 14a also illustrates the distribution loss of 
block attributes after Segment 0 is lost during transmission. In 
particular, Loss Pattern 1415 shows the DR, Motion Flag, and MIN 
values loss across six segments after a subsequent deshuffling is 
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applied to the received data that was initially shuffled using Modular 
Shuffles 1410, 1420, and 1430. As illustrated in Figure 14a, the block 
attribute loss is distributed periodically across Segments 0, 2, and 4 
while Segments 1, 3, and 5 have no block attribute loss. Additionally, 
Spatial Loss Pattern 1417 illustrates the deshuffled spatial distribution 
of damaged FL-data after Segment 0 is lost during transmission. 
Spatial Loss Pattern 1417 shows the DR, Motion Flag, and MIN value 
loss after a subsequent deshuffling is applied to the received data. In 
Spatial Loss Pattern 1417, a damaged block is surrounded by 
undamaged blocks and damaged block attributes can be recovered with 
surrounding undamaged blocks. 

Figure 14 and Figure 14a illustrate a modular three shuffling 
pattern and the distribution loss of block attributes after a segment is 
lost during transmission. In alternative embodiments, the count 
variables or the number of segments are varied to alternate the 
distribution of lost block attributes. Figure 14b illustrates Modular 
Shuffle Result 1421 and Loss Pattern 1420. Similarly, Figure 14c 
illustrates Modular Shuffle Result 1426 and Loss Pattern 1425. Both 
Loss Pattern 1420 and Loss Pattern 1425 illustrate the distribution loss 
of block attributes across six segments, as opposed to three segments as 
previously described. 

It is contemplated that in alternate embodiments various 
combinations of block attributes will be distributed to perform the 
shuffling process. 

Inter Segment VL-Data Shuffling 
In the inter segment VL-data shuffling process, bits between a 
predetermined number of segments, for example, 6 segments, are 
arranged to ensure a spatially separated and periodic VL-data loss 
during an up to 1/6 packet transmission loss. Figure 15 and 16 
illustrate one embodiment of the inter segment VL-data shuffling 
process. 

In the present embodiment, a transmission rate approaching 30 
Mbps is desired. Accordingly, the desired transmission rate results in 
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31,152 bits available for the VL-data in each of the 60 buffers. The 
remaining space is used by FL-data for the eighty eight blocks included 
in a buffer. Figure 15 includes the VL-data buffer organization within 
a frame set for a transmission rate approaching 30 Mbps. As 
previously described, partial buffering is used to maximize the usage 
of available VL-data space within each buffer, and the unused VL-data 
space is filled with a post-amble. 

Figure 16 illustrates one embodiment of the shuffling process to 
ensure a spatially separated and periodic VL-data loss. The first row 
illustrates the VL-data from the 60 buffers in Figure 15 rearranged into 
a concatenated stream of 1,869,120 bits. The second row illustrates the 
collection of every sixth bit into a new stream of bits. Thus, when the 
decoder subsequently reverses the process, a burst loss of up to 1/6 of 
the data transmitted is transformed into a periodic loss where at least 5 
undamaged bits separate every set of two damaged bits. 

The third row illustrates grouping every 10 bits of Stream 2 into 
a new stream of bits, Stream 3. The boundary of a grouping is also 
defined by the number of bits in a segment. Grouping of Stream 2 for 
every tenth bit ensures that a 1/60 data loss results in fifty-nine 
undamaged bits between every set of two damaged bits. This provides 
for a spatially separated and periodic VL-data loss in the event that 88 
consecutive packets of data are lost. 

The fourth row illustrates grouping every 11 bits of Stream 3 
into Stream 4. The boundary of a grouping is also defined by the 
number of bits in a segment. Grouping of Stream 3 for every eleventh 
bit ensures that 1 / 660 data loss results in 659 undamaged bits between 
to damaged bits, resulting in a spatially separated and periodic VL-data 
loss during a transmission loss of 8 consecutive packets. 

Each group of 31,152 bits within Stream 4 is consecutively re- 
stored in Buffers 0 - 59, with the first group of bits stored in Buffer 0 
and the last group of bits stored in Buffer 59. 

It will be appreciated by one skilled in the art that the grouping 
requirements of Figure 16 are variable to ensure a spatially separated 
and periodic VL-data loss tolerance up to a 1/n transmission loss. 
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Transmission 

The previously described shuffling process creates buffers with 
intermixed FL-data and VL-data. For one embodiment, packets are 
generated from each buffer, according to packet structure 200, and 
transmitted across Transmission media 135. 

Data Recovery 

As noted earlier, the innovative method for encoding the 
bitstream of data enables robust recovery of data that typically occurs 
due to lost packets of data. The general overview of the decoding 
process has been shown in Figure 4. 

Referring to Figure 4, the data received in packets is processed 
through the multiple level deshuffling process, steps 425, 430, 435, and 
440 wherein different levels or portions of the bitstream received via 
packets are deshuffled to retrieve data. ADRC decoding is then 
applied to the data, step 445, in accordance with the teaching known in 
the art (e.g., Kondo, Fujimori, Nakaya, "Adaptive Dynamic Coding 
Scheme for Future HDTV Digital VTR", Fourth International 
Workshop on HDTV and Beyond, Sept. 4 - 6, 1991, Turin, Italy). 

Intra frame set block deshuffling is then performed and block- 
to-image mapping is subsequently executed, steps 450, 455. Steps 425, 
430, 435, 440, 445, 450, and 455 are inverse processes of the earlier 
process steps performed to encode the data and will not be discussed in 
detail herein. However, it should be noted that in one embodiment, 
deshuffling levels represented by steps 425, 430 and 440 are data 
independent. For example, the deshuffling process performed is 
predetermined or specified by an address mapping or table lookup. 
Since deshuffling steps 425, 430 and 440 are independent of data 
contents, data loss due to, for example, packet loss, does not prevent 
the deshuffling steps from being performed. Similarly, steps 450 and 
455 are data independent. The intra group VL-data deshuffling 
process, however, is dependent on the contents of data. More 
particularly, the intra group VL-data deshuffling process is used to 
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determine the quantization codes for the blocks of the groups. Thus, 
at step 435, if packets are lost, the affected groups cannot be processed. 

After execution of the deshuffling, decoding and mapping (steps 
425, 430, 435, 440, 445, 450 and 455), a recovery process is performed to 
recover the Qbit and Motion Flag values that were located in lost 
packets. The Qbit value is lost typically due to DR loss (due to lost 
packets). When the Qbit or Motion Flag value is unknown, the Q 
code bits of a pixel cannot be determined from the data bitstream. If a 
Qbit or Motion Flag value is improperly determined then this error 
will propagate as the starting point of subsequent blocks in that data in 
the buffer will be incorrectly identified. 

Figure 17 describes the general process for recovering the Qbit 
and Motion Flag values. This particular embodiment describes the 
process using multiple blocks of data to recover the Qbit and Motion 
Flag values; however, it is contemplated that the particular number of 
blocks is not limited by the discussion herein and could be one or 
more blocks. Referring to Figure 17, based on the detection of an error 
in the bitstream, step 1705, candidate decodings based on specified 
parameters are generated for the three blocks examined. At step 1715, 
each candidate decoding is scored on the likelihood that it is an 
accurate decoding and at step 1720, the candidate decoding with the 
best score is used, the decoding identifying the Qbit and Motion Flag 
values which enable the subsequent decoding of pixels of the affected 
blocks. 

Referring back to the decoding process of Figure 4, once the best 
decoding is selected, any DR or MIN values that were lost due to lost 
packets are recovered, step 465. A variety of recovery processes known 
to one skilled in the art can be applied to recover DR and MIN, 
including least squares or the averaging of values of adjacent blocks. 
For one example, see, Kondo, Fujimori, Nakaya, "Adaptive Dynamic 
Coding Scheme for Future HDTV Digital VTR", Fourth International 
Workshop on HDTV and Beyond, Sept. 4 - 6, 1991 Turin, Italy. In the 
present embodiment, the innovative image-to-block mapping 
processes and data structures created therefrom increase the number 
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of neighboring blocks, thus providing additional data and facilitating 
more accurate DR or MIN recovery. In particular, in one 
embodiment, DR and MIN are recovered as follows: 



m-X [ yi -MINy qi 
DR = f ^ , 



where DR' corresponds to the recovered DR, qj is the i-th value in an 
ADRC block and q s e {0,1,...2Q-1}; m = 2<M for Edge-matching ADRC 
and m = 2<2 for Non-edge-matching ADRC; yi is a decoded value of an 
adjacent block pixel; and Q is the Qbit value; and 



MIN - 



N 

where MIN' corresponds to the recovered MIN and N is the number 
of terms used in the summation (e.g., N=32 when i = 0-31). In another 
embodiment, if DR and MIN of the same block are damaged at the 
same time, DR and MIN are recovered according to the following 
equations: 



m 

DR = — 



MIN 

N 



At step 470, ADRC decoding is applied to those blocks not 
previously decoded prior to Qbit and Motion Flag recovery and a pixel 
recovery process is executed, step 475, to recover any erroneous pixel 
data that may have occurred due to lost packets or random errors. In 
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addition a 3:1:0 -> 4:2:2 back conversion is performed, step 480, to place 
the image in the desired format for display. 

Figure 18 illustrates one particular embodiment of the Qbit and 
Motion Flag recovery process of the decoding process of the present 
invention. In this particular embodiment, the inputs to the process 
are adjacent block information, and the block attributes and pixel data 
for the three blocks to be processed. Error flags indicating the location 
of the lost data are also input. The error flags can be generated in a 
variety of ways known to one skilled in the art and will not be 
discussed further herein except to say that the flags indicate which bits 
were transmitted by damaged or lost packets. 

At step 1805, the candidate decodings are generated. The 
candidate decodings can be generated a variety of ways. For example, 
although the processing burden would be quite significant, the 
candidate decodings can include all possible decodings. Alternately, 
the candidate decodings can be generated based on pre-specifted 
parameters to narrow the number of candidate decodings to be 
evaluated. 

In the present embodiment, the candidate decodings are 
determined based on the possible key values used to randomize a 
bitstream of the intra group VL-data shuffling process earlier 
described. In addition, it should be noted that candidate decodings are 
further limited by the length of the bits remaining to be decoded and 
knowledge of how many blocks remain. For example, as will be 
discussed, if processing the last block typically the decoding length of 
that block is known. 

Continuing with the present example, Figure 19 illustrates 
possible cases for the present embodiment where the value x indicates 
an unknown value (which may be due to packet loss). This is further 
explained by example. m\ is defined as the Motion Flag of the i-th 
block, qi is the number of the quantization bits of the i-th block, n; is 
the number of possible candidates of the i-th block and dj is the value 
of a key element of the i-th block described previously in intra group 
VL-data shuffling. The i-th block is defined within each group. In this 
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example, the number of blocks within each group is three. A key for 
the three block group is generated as, do + 10-di + 100 -d2. Assuming 
that in the first block the Motion Flag is unknown and the number of 
quantization bits is 2, mo equals x and qo equals 2. Following the 
equation described above to generate the key element, dj = 5 m j 
+ qi, the set of possible digits for do consists of {2 and 7}. Thus, the 
number of possible values (no) is 2. Assuming the second block to 
have a Motion Flag value of 1 and one quantization bit, and the value 
for di is 5-1+1 = 6 and ni = 1. The third block has a Motion Flag value 
of 1 and an unknown number of quantization bits. Thus, the digit d2 
includes a set consisting of {6, 7, 8, 9} and n2 = 4. Thus, the number of 
possible candidates of this group, M, is 2 • 1 • 4 = 8, and the keys used to 
generate the candidate decodings are the variations of 662, 667, 762, 
767, 862, 867, 962, 967. This process is preferably used for each group 
which was affected by data loss. 

Referring back to Figure 17, at step 1715, once the data has been 
decoded in accordance with the key data, the candidate decodings 
generated are evaluated or scored on the likelihood that it is a correct 
decoding of the data. At step 1720, the candidate decoding with the 
best score is selected to be used. 

A variety of techniques can be used to score the candidate 
decodings. For example, the score may be derived from an analysis of 
how pixels of blocks of a particular candidate decoding fit in with 
other pixels of the image. Preferably the score is derived based upon a 
criteria indicative of error, such as a square error and correlation. For 
example, with respect to correlation, it is a fairly safe assumption that 
the adjacent pixels will be somewhat closely correlated. Thus, a 
significant or a lack of correlation is indicative that the candidate 
decoding is or is not the correct decoding. 

As is shown in Figure 18, four different criteria are analyzed to 
select the best candidate decoding. However, it is contemplated that 
one, two, three or more different criteria can be analyzed to select the 
best candidate decoding. 
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Referring to Figure 18, the present embodiment utilizes four 
subscoring criteria which are subsequently combined into a final score. 
In particular, in step 1815, the square error measure is generated, step 
1820, horizontal correlation is determined, step 1825, vertical 
correlation is determined, and at step 1830 temporal activity is 
measured (each an M-by-2-N matrix in accordance with M candidates, 
N blocks and 2 frames/block of data). Although horizontal and 
vertical correlation is used, it should be recognized that a variety of 
correlation measurements, including diagonal correlation, can be 
examined. At steps 1835, 1840, 1845, 1850, a confidence measure is 
generated for each criterion to normalize the measurements 
generated, and at steps 1855, 1860, 1865 and 1870, a probability function 
for each of the different criteria is generated. These probability 
functions are then combined, for example, by multiplying the 
probability values to generate a score, for example, the likelihood 
function shown in Figure 18, step 1875. The score for the candidate 
decoding is subsequently compared against all candidate decoding 
scores to determine the likely candidate. 

It should be recognized that a variety of techniques can be used 
to evaluate the candidate decodings and generate the "scorings" for 
each candidate. For example, confidence measures are one way of 
normalizing the criteria. Furthermore, a variety of confidence 
measures, besides the ones described below, can be used. Similarly, 
multiplying the probability values based on each criterion to generate 
a total likelihood function is just one way of combining the variety of 
criteria examined. 

The encoding processes facilitate the determination of the best 
candidate decoding because typically the candidate decodings which 
are not the likely candidate, will have *a relatively poor score, while 
decodings that are quite likely candidates will have a significantly 
better score. In particular, the Q code randomization process described 
previously in the intra group VL-data shuffling process assists in this 
regard. 
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Figures 20a 20b, 20c and 20d provide illustrations of the 
different measurements performed at steps 1815, 1820, 1825 and 1830 of 
Figure 18 to generate the scoring and total score for a particular 
candidate decoding. Figure 20a illustrates the square error to evaluate 
a candidate decoded pixel x { as compared to its decoded neighbors y^, 
wherein the suffix "i,j" is corresponding to the neighboring address of 
T\ It is preferred that some of the largest terms are removed to 
remove any influences due to spikes, that is the terms that arise due to 
legitimate edges in the image. Preferably, the three largest terms of (x r 
y { j) 2 are discarded to remove spikes that may occurred. Figure 20b 
illustrates the temporal activity criteria. This is applicable only when 
it is or is assumed to be a motion block. The temporal activity criteria 
assumes that the better the candidate decoding, the smaller the 
differences between blocks. Thus the worse the candidate decoding, 
the larger the differences between blocks. Spatial correlation assumes 
that the more likely candidate decodings will result in heavy 
correlations as real images tend to change in a slow consistent way. 
The horizontal correlation process illustrated in Figure 20c and 
vertical correlation process illustrated by Figure 20d utilize that 
assumption. 

The confidence measures, steps 1835, 1840, 1845, and 1850 of 
Figure 18, provide a process for normalizing the criteria determined in 
the previous steps (steps 1815, 1820, 1825 and 1830). In one 
embodiment, for example, the confidence measure for the square 
error takes values from the interval [0,1], and confidence is equal to 0 
if the errors are equal and equal to 1 if one error is 0. Other measures 
or methods to normalize are also contemplated. 

Similarly, the confidence measure for the spatial correlation is: 

maximum(Y,0) - maximum(X,0) 
where Y is the best correlation value and X is the correlation for the 
current candidate decoding. The temporal activity confidence 
measure is determined according to the following equation: 

conf = (a-b) / (a+b) 
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where a = max (X, M_TH) and b = max (Y,M_TH) where M_TH is the 
motion threshold for the candidate block and Y is the best 
measurement, that is the smallest temporal activity, and X equals the 
current candidate measurement of temporal activity. 

At steps 1855, I860, 1865 and 1870, Figure 18, the probability 
function is generated for each of the different criteria. A variety of 
methods can be used to generate the probability measure. For 
example, a score can be prescribed to a confidence measure. If the 
confidence measure is greater than a predetermined value, e.g., 0.8, 
the base score is decreased by 10; if between 0.5 and 0.8, the base score 
decreased by 5. Figure 21 illustrates one embodiment in which a table 
used to generate the probability function for the square error 
measurement criteria. The table includes empirically determined 
based data arbitrarily binned containing confidence and square error 
measures and known candidate decodings. More particularly, the 
table can be generated by using undamaged data and assuming that the 
DR was corrupted or lost. Keys and confidence measures for correct 
and incorrect decodings are then generated. The table reflects the 
probability ratio of correct to incorrect decodings. Using this table, for 
a particular squared error value (row) and confidence value (column), 
the probability can be determined. For example, it can therefore be 
seen that for a variety of square error measures at a confidence 
measure of zero, there is approximately a 40% to 50% probability that 
the candidate is correct. If the confidence is not 0, but small, the 
probability drops significantly. Similar probability tables are generated 
for the correlation and temporal measurements based on 
corresponding empirically determined criteria measurements and 
confidence measurements. 

The probabilities generated are considered data to generate 
"scores" in the present embodiment and as noted earlier, other 
techniques to score candidate decodings may be used. At step 1875, the 
different probabilities are combined into a likelihood function Lj = 
j -Pi,j, where j is a multiplication function of probability functions 
Pi,j, and Pi,j,is the probability function for candidate i, block j. The 
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candidate is therefore selected as the one that maximizes the function 
U 

Referring back to Figure 18, it may be necessary to recover 
certain block attributes that were transmitted in lost packets. 
Therefore, at step 1810, DR and MIN values are recovered where 
necessary. A variety of techniques, from default values, averaging, 
squared error functions to more sophisticated techniques, including 
those discussed in Kondo, Fujimori and Nakaya, "Adaptive Dynamic 
Range Coding Scheme for Future HDTV Digital VTR", and Kondo, 
Fujimori, Nakaya and Uchida, "A New Concealment Method for 
Digital VCRs", IEEE Visual Signal Processing and Communications, 
September 20-22, 1993, Melbourne Australia, may be used. The 
recovered values are utilized to generate the candidate decodings as 
discussed above. 

Alternately, the DR and MIN values are determined during the Qbit 
determination process. This is illustrated in Figure 22. In particular, 
as noted above, in the present embodiment, the Motion Flag and 
number of quantization bits are used in the encoding process and later 
used during the recovery process to narrow the number of possible 
candidate decodings. As noted earlier, other information can also be 
used. Thus the value of DR and /or value of MIN may also be used to 
encode the data. Alternately, a portion of bits of DR are used for 
encoding (e.g., the two least significant bits of DR). Although the DR 
data is encoded, the number of possible candidate decodings is 
increased significantly as variables are added. Referring to Figure 22, 
KM candidate decodings are therefore generated, where K is the 
number of candidate values for the unknown data, e.g. K=4 if two bits 
of the sum of DRi, DR 2 and DR 3 is encoded (DRi, DR 2 and DR3 
represent the DR values of the blocks of the group). The DR and MIN 
are therefore recovered using the auxiliary information provided, e.g., 
the encoded two bits of the sum of DRi, DR2 and DR3. This improves 
the process of candidate selection at the cost of additional overhead to 
examine the larger number of candidate decodings. 
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It should be noted that generally, the more neighboring blocks 
that are decoded, the better the Qbit and Motion Flag recovery process. 
Furthermore, in some embodiments the process is applied to each 
subsequent block of a buffer; if ail or some of the FL-data is available, 
the number of candidate decodings can be reduced, possibly to one 
candidate decoding given all the FL-data for a block is available. 
However, it is desirable that the Qbit and Motion Flag recovery process 
be avoided altogether as the process is a relatively time consuming 
one. Furthermore, it is desirable to use as much information as 
possible to perform Qbit and Motion Flag recovery. In one 
embodiment, blocks are processed from the beginning of a buffer until 
a block with lost Qbit/ Motion Flag information is reached. This is 
defined as forward Qbit and Motion Flag recovery. In another 
embodiment, the end of the buffer is referenced to determine the 
location of the end of the last block of the buffer and the data is 
recovered from the end of the buffer until a block with lost 
Qbit/Motion Flag data is reached. This is defined as backward Qbit and 
Motion Flag recovery. 

As noted earlier, the blocks are variable in length, due the 
length of the VL-data; therefore there is a need to determine the 
number of bits forming the VL-data of a block so that the position of 
subsequent blocks in the buffer can be accurately located. During the 
encoding process, a post-amble of a predetermined and preferably 
easily recognizable pattern is placed in the buffer to fill the unused bit 
locations. During the decoding process, the post-amble will be located 
between the block and the end of the buffer. As the pattern is one that 
is easily recognizable, review of patterns of bits enables the system to J 
locate the beginning of the post-amble and therefore the end of the last 
block in the buffer. This information can be used in two ways. If the 
last block contains damaged Qbit/Motion Flag data and the beginning 
of the last block is known (e.g., the preceding blocks have been 
successfully decoded), the difference between the end of the 
immediate preceding block and the beginning of the post-amble 
corresponds to the length of the block. This information can be used 
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to calculate the Qbit and /or Motion Flag of the block. The starting 
location of the post-amble can also be used to perform Qbit and 
Motion Flag recovery starting at the last block and proceeding towards 
the beginning of the buffer. Thus, the Qbit and Motion Flag recovery- 
process can be implemented bidirectionally. 

Figure 23 illustrates the use of a post-amble in the bidirectional 
Qbit and Motion Flag recovery process. Referring to Figure 23, the 
buffer 2300 includes FL-data 2303 for the N groups of blocks of VL-data. 
Each group consists of a plurality of blocks (e.g., 3 blocks). In the 
present example, the first two groups 2305, 2310 are decoded and the 
third group 215 cannot immediately be decoded due to damaged 
DR/Motion Flag data. At this point, the Qbit/ Motion Flag recovery 
process is required in order to recover the damaged data. Rather than 
continue processing groups in the forward direction, the process refers 
to the end of the buffer, determined by looking for the post-amble 
pattern 220. The beginning of the post-amble and therefore the end of 
the last group of blocks are determined. As the DR/Motion Flag data 
is indicative of the length of the VL-data, the beginning of the VL data 
of the last block, and therefore the end of the immediate preceding 
block, is determined. Therefore, the blocks can be decoded , e.g., 
blocks 225, 230, 235 until a block 240 with damaged data is reached. 
The damaged 215, 240 and obstructed blocks 250 are then recovered, 
preferably using the Qbit/Motion Flag recovery process described 
above. 

It should be noted that the bidirectional process is not limited to 
a sequence of forward and reverse processing; processing can occur in 
either or both directions. Furthermore, in some embodiments, it may 
be desirable to perform such processing in parallel to improve 
efficiency. Finally, it is contemplated that undamaged obstructed 
blocks may be recovered by directly accessing the Qbit/Motion Flag 
information without executing the Qbit/Motion Flag recovery process 
described above. 

As noted earlier, a variety of scoring techniques may be used to 
determine the best candidate decoding to select as the decoding. In an 
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alternate embodiment, the smoothness of the image using each 
candidate decoding is evaluated. In one embodiment, the Laplacian 
measurement is performed. The Laplacian measurement measures a 
second-order image surface property, e.g., surface curvature. For a 
linear image surface, i.e., smooth surface, the Laplacian measurement 
will result in a value that is approximately zero. 

The process will be explained with reference to Figures 24a, 24b, 
and 24c. Figure 24a illustrates one embodiment of the Laplacian 
kernel. It is contemplated that other embodiments may also be used. 
The kernel "L" represents a 3x3 region. To measure smoothness of 
the image region, 3x3 subregions of the image (Figure 24b) are 
convolved with the kernel and the convolved values are averaged. 
The size of the region and subregion (and therefore kernel size) can be 
varied according to application. 

One embodiment of the process is described with reference to 
Figure 24c. This embodiment utilizes a kernel and subregion size of 
3x3 and a region size of 8x8, the individual elements identified by 
indices i,j. At step 2460, the candidate decoded values x[i][j] are 
normalized. For example, the values can be normalized according to 
the following equation: 

x[i][j] = -, x[iM o<a,j<8 



^Z(x[i][j]-X mMn ) : 



where, = '•> ^ , 0</,;<8 

At step 2465, the normalized values are used to compute a block 
Laplacian value Lx indicative of smoothness according to the 
following: 

l[i][j]= E I L[m][n] x[i + m][j+n], 0<a,j<8 
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— iJ 



x 64 

The closer the block Laplacian value is to zero, the smoother the 
image portion. Thus a score can be measured based upon the block 
Laplacian value, and the decoding with the least Laplacian value is the 
correct one. 

The Laplacian evaluation can also be achieved using candidate 
encoded values q[i][j] . The basic process is the same as the candidate 
decoded value case of Figure 24c. This embodiment utilizes a kernel 
and subregion size of 3x3 and a region size 8x8, the individual 
elements identifies by the indices i,j. At step 2460, the candidate 
encoded values q[i][j] are normalized. For example, the values can be 
normalized according to the following equation: 



MM 

where, Q mean = -^rj 

At step 2465, the normalized values are used to compute the block 
Laplacian value L q indicative of smoothness according to the 
following equation: 

'['!/■]= i X4>"M<7'['- + H[; + 4 l*i.y<7 

m=-l n=— I 



L = 

q 36 



The closer the block Laplacian value is to zero, the smoother the 
image portion. Thus a score can be measured based upon the block 
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Laplacian value and the candidate with the smallest Laplacian value is 
the correct one. 

Other variations are also contemplated. In alternative 
embodiments, higher order image surface properties can be used as a 
smoothness measure. In those cases, higher order kernels would be 
used. For example, a fourth order block Laplacian measurement may 
be performed using a fourth order kernel. Such a fourth order kernel 
can be realized using two second order Laplacian computations in 
cascade. 

It is further contemplated that the evaluation process is 
dependent upon whether the image has an activity or motion larger 
than a predetermined level. If the image portion is evaluated to have 
larger motion than a predetermined level, then it is preferable to 
perform the measurements on a field basis as opposed to on a frame 
basis. This is explained with reference to Figure 25. Figure 25 explains 
the process using smoothness measures; however, it is contemplated 
that this process can be implemented using a variety of types of 
measures. 

Frame 2505 of an image region is composed of field 0 and field 
1. If motion is not detected, step 2510, the smoothness measurement 
is computed by computing the block Laplacian value for the block 
within each frame, step 2515. If larger motion than a predetermined 
level is detected, block Laplacian measurements are performed on 
each field, steps 2520, 2525 and the two measurements are combined, 
step 2530, , e.g., averaged to generate the smoothness measurement. 

Motion can be detected /measured a variety of ways. In one 
embodiment, the extent of change between fields is evaluated and 
motion is detected if it exceeds a predetermined threshold. 

Motion detection and the use of frame information and field 
information to generate recovered values (typically to replace lost or 
damaged values) can be applied to any portion of the process that 
requires a recovered value to be generated. For example, motion 
detection and the selective use of frame information and field 
information to generate recovered values can be applied to DR/MIN 
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recovery, pixel recovery as well as Qbit and -Motion Flag recovery 
processes. Thus, based on the level of motion detected, the recovery 
process will utilize existing information on a field basis or frame basis. 
Furthermore, this process can be combined with the application of 
weighting values that are selected based upon levels of correlation in 
particular directions (e.g., horizontal or vertical). 

In another embodiment of the Qbit and Motion Flag recovery 
process, candidate decodings are evaluated based upon intra block and 
inter block measurements. In the following discussion, the term 
"block" refers to a portion of a frame or field. The intra block 
measurement evaluates the candidate decoded image portion, e.g., the 
smoothness of the image portion. The inter block measurement 
measures how well the candidate decoding fits with the neighboring 
image portions. Figures 26a and 26b illustrate the combined inter 
block and intra block evaluation. In particular, Figure 26a shows an 
acceptable candidate decoding as both the inter block and intra block 
measurements are good, whereas in Figure 26b the inter block 
measurement is poor, even though the intra block measurement is 
quite good. 

Examples of intra block measurements include the smoothness 
measurement described above. Examples of inter block 
measurements include the square error measurements described 
earlier. An alternative inter block measurement is the ratio of 
compatible boundary pixels and the total number of boundary pixels at 
the candidate ADRC block. 

An example of an inter block and intra block evaluation of an 
8x8 block that is ADRC encoded will be explained with respect to 
Figures 26c, 26d and 26e. Figure 26d illustrates an image portion 
(block) of data of a encoded values 2650 consisting of q values from 
which candidate decoded values x are generated and neighboring 
decoded data 2655 consisting of y values. As set forth in the flow chart 
of Figure 26c, at step 2605, the intra block measure is computed to 
generate a measure, e.g., block Laplacian L x . At step 2610, the inter 
block measure S x is computed to generate a measure of compatibility 
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between adjacent blocks. At step 2615, the combined measure M x is 
generated. The combined measure provides the information used to 
select a candidate decoding. 

In the present embodiment, S x is computed as the number of 
neighboring data that lies in a valid range for each boundary pixel of 
candidate decoding (see Figure 26e). Figure 26e is a chart illustrating a 
valid range for one embodiment which shows a valid range of each 
observed quantized value q*. Thus Lq DR < Uq, where Lq, Uq 
respectively represent the lower and upper bounds of DR 
corresponding to the number of quantization bits = Q. Preferably S x is 
normalized according to the following: S x = S x /number of boundary 
pixels. 

In the present embodiment the combined measure M x is 
computed according to the following equation: M x = S x + (1-L X ). 
Alternatively, the combined measure may be weighted such that the 
following equation would be used: Mx =w-S x + (1 - w)(l-L x ), where w 
is the weighting value, typically an empirically determined weighting 
value. 

Other embodiments for determining DR and MIN values that 
have been lost/damaged are also contemplated. For example, the 
earlier described equations can be modified to recover DR and MIN 
values with higher accuracy. In an alternate embodiment, a median 
technique is applied. In one embodiment of the median technique, 
the value of MIN is recovered as the median of all MINi values 
computed as: 

MINi = yi-qi • s 

where qi represents the encoded pixel value and yi represents the 
decoded pixel neighboring qi. For edge-matching ADRC, s = DR/(2Q - 
1). For non-edge-matching ADRC, s = DR/2Q, where Q represents the 
number of quantization bits per pixel (Qbit value). 

The values used may be temporally proximate or spatially 
proximate. The values of yi may be the decoded value of the 
neighboring pixel in an adjacent frame/field or the same field. The 
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values of yi may be the decoded value of the pixel from the same 
location as qj in an adjacent frame/field or the same field. 

In addition, any DR and /or MIN recovery technique may be 
combined with a clipping process to improve recovery accuracy and 
prevent data overflow during the recovery process. The clipping 
process restricts the recovered data to a predetermined range of values; 
thus those values outside the range are clipped to the closest range 
bound. In one embodiment, the clipping process restricts values in 
the range [Lq, Uq], where Lq, Uq respectively represent the lower and 
upper bounds of the range of pixel values represented by the number 
of quantization bits = Q. quantization bits, and further restricts values 
to: MIN + DR Num, where Num represents the maximum pixel 
value; in the present embodiment, Num is 255. In the present 
embodiment, where applicable, Uq +1 = Lq+i 

Combining the criteria into a single equation results for an 
unbounded recovered value (val*) for the DR, the final clipped 
recovered value (val) is obtained from the following equation: 

val = max(min(val, min(UQ,255-MIN)),LQ) 
where min and max respectively represent minimum and maximum 
functions. 

In an alternate embodiment, the boundary pixels yi used to 
generate an recovered DR and/or MIN can be filtered to only use 
those that appear to correlate best, thereby better recovering DR and 
MIN. Those boundary pixels not meeting the criteria are not used. In 
one embodiment, a boundary pixel yi is considered valid for DR 
calculations if there exists a value of DR such that Lq DR<Uq and an 
original pixel yi would have been encoded as qi. Thus, a pixel is valid 
if the following equations are satisfied: 
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(y t -MIN)m ^ 
max(^ -0.5,0) - Q 

( yi -MIN)m 
min(^ J r — 0.5,m) 

where m represents the maximum quantization level = 2Q-1. A DR 
recovered value (val*) can then be computed according to the 
following equation: 



vaf = 



The value can then be clipped into the valid range. Thus this 
process forces the DR recovered value into the interior of the valid 
region as defined by the threshold table, reducing the accuracy for 
points whose true DR lies near the threshold table boundary. 

It has been noted that due to quantization noise, the DR of 
stationary ADRC blocks varies slightly from frame to frame. If this 
variance crosses an ADRC encoding boundary, and if the DR is 
recovered on several consecutive frames, then the DR recovered 
value with valid pixel selection tends to overshoot at each crossing, 
resulting in a noticeable blinking effect in the display. In an attempt to 
reduce the occurrence of this effect, in one embodiment, the valid 
pixel selection process is modified to relax the upper and lower 
bounds, allowing border pixels that encroach into the neighboring 
valid region. By including points just outside the boundary, it is more 
likely that the recovered value will take on a value near that of the 
upper or lower bound. The relaxed bounds L'q and U'q are computed 
by means of a relaxation constant r. In one embodiment, r is set to a 
value of .5. Other values can be used: 
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L'q = tLq-i + (1-r) Lq 
U'q = (l-r)U Q + rU Q+ i 

The discussion above sets forth a number of ways to recover DR 
and MIN when the values have been damaged or lost. Further 
enhancements can be realized by examining the correlation between 
data temporally and /or spatially, and weighting corresponding 
calculated recovered values accordingly. More particularly, if there is 
a large correlation in a particular direction or across time, e.g., 
horizontal correlation, there is a strong likelihood that the image 
features continue smoothly in that direction that has a large 
correlation and therefore an recovered value using highly correlated 
data typically generates a better estimate. To take advantage of this, 
boundary data is broken down into corresponding directions (e.g., 
vertical, horizontal, field- to-field) and weighted according to the 
correlation measurement to generate a final recovered value. 

One embodiment of the process is described with reference to 
Figure 27a. At step 2710, a recovered value of the DR or MIN value to 
be recovered is generated in one direction and at step 2715, a recovered 
value is generated in another direction. For example, if the process is 
spatially adaptive, then boundary pixels along horizontal borders are 
used to generate a first recovered value, "hest", and boundary pixels 
along vertical borders are used to generated a second recovered value, 
"vest". Alternately, if the process is temporally adaptive, then 
boundary pixels between adjacent fields are used to generate a first 
recovered value and boundary pixels between adjacent frames are 
used to generate a second recovered value. 

At step 2720, the recovered values are weighted according to 
correlation calculations indicative of the level of correlation in each 
direction. The weighted first and second recovered values are 
combined to generate a combined recovered value, step 2725. It 
should be noted that the process is not limited to generated weighted 
recovered values in only two directions; it is readily apparent that the 
number of recovered values that are weighted and combined can be 
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varied according to application. A variety of known techniques can be 
used to generate a correlation value indicative of the level of 
correlation in a particular direction. Furthermore, a variety of criteria 
can be used to select the weighting factor in view of the levels of 
correlation. Typically, if one correlation is much larger than the other, 
the combined recovered value should be based primarily on the 
corresponding recovered value. In one embodiment, the combined 
recovered value is computed as follows: 



vaf-i 



ahest + (l- a)vest : hc>vc 
(1 - a)hest + ccvest : he < vc 



where he represents the horizontal correlation, vc represents the 
vertical correlation, hest represents a DR recovered value based only 
on left and right boundary information, and vest represents a DR 
recovered value based only on top and bottom boundary information, 
and a represents, the weighting value. The weighting value can be 
determined a variety of ways. Figure 27b illustrates one embodiment 
for determining weighting values as a function of the difference 
between the horizontal correlation and vertical correlation. More 
particularly, a was chosen to be: 



a(\hc - vc\) = 



0.5 + 0. 25 -£" 8(0 - 35 "" ,c " vc,) : \hc-vc\< 0.35 
l-0.25.e- 8(,Ac - wH>35) : |*c- vc| > 0.35 



As noted above, the adaptive correlation process is applicable to 
both DR and MIN recovery. It is preferred, however, that the MIN 
recovery is clipped to insure that MIN + DR 255, therefore the 
function val = max(min(val', 255 - MIN), 0) can be used. 
Furthermore, as noted above, the temporal correlation can be 
evaluated and used to weight recovered values. In addition, a 
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combination of temporal and spatial correlation can be performed. 
For example, one recovered value is generated between fields as a 
temporal recovered value. Another recovered value is generated 
within one field as a spatial recovered value. The final recovered 
value is computed as the combination value with a combination of 
temporal and spatial correlation. The correlation combination can be 
replaced with a motion quantity. Other variations are also 
contemplated. The method can also be applied to audio data. 

In an alternate embodiment, a low complexity modification to 
the least squares technique is used. Using this embodiment the 
blinking experienced due to recovered DR values is reduced. For 
purposes of the following discussion, QV represents a list of encoded 
values from the image section or ADRC block whose DR is being 
recovered having a set of points qi and Y is a list of decoded values 
taken from the vertical or horizontal neighbors of the points in QV, 
where yi represents a vertical or horizontal neighbor of qi. As each 
point qi may have up to four decoded neighbors, one pixel or point 
may give rise to as many as four (qi, y\) pairings. The unconstrained 
least squares estimate of DR (DR u i s ) is thus: 

2 fl -£(y,-AfflV)-(0.5 + ft ) 
iDR)u "= ' X(0-5 + ft ) a 

where Q is the number of quantization bits, MIN is the minimum 
value transmitted as a block attribute. The above equation assumes 
non-edge-matching ADRC; for edge-matching ADRC, 2Q is replaced 
with 2Q-1 and (0.5 + qi) is replaced with qi. 

The unconstrained least squares estimate is preferably clipped to 
assure consistency with the threshold table and the equation MEN + 
DR 255 which is enforced during encoding (Typically, for non-edge- 
matching ADRC, permissible DR values are in the range of 1-256). 
Thus, the least squares estimate is clipped (DRi sc ) by: 
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(DR)i S c = max(min(UB,DR u is),LB) 

where UB represents the upper bound and LB represents the lower 
bound and min and max respectively represent minimum and 
maximum functions. 

In an alternate embodiment, the estimation can be enhanced by 
selecting the pixels that are more suitable for DR estimation to 
calculate the estimate of DR. For example, flat regions in an image 
provide pixels which are more suitable for DR estimation than those 
regions in which high activity occurs. In particular, a sharp edge in 
the edge may decrease the accuracy of the estimate. The following 
embodiment provides a computationally light method for selecting 
the pixels to use to calculate an estimate of DR. 

In one embodiment, the least squares estimate (DRi se )/ e.g., 
DR u i s or DRi sc . is computed. Using this estimate, the list of encoded 
values QV is transformed into candidate decoded values X, where xi 
are members of X derived from qi. The Xi value is a recovered 
decoded value formed using the first estimate of DR. The X* value is 
defined according to the following equation: 

Edge-matching ADRC: x t . = MIN + (o. 5 + 



Non-edge-matching ADRC: jc. = MIN + 




Assuming DRi se is a reasonable estimate of the true DR, then 
anywhere that Xi is relatively close to yi, may be judged to be a low 
activity area and thus a desirable matching. New X and Y lists may 
then be formed by considering only the matches where xi and yi are 
close and the least squares estimate recomputed to generate an 
updated estimate. 
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The criteria for determining what is considered "close" can be 
determined a number of ways. In one embodiment, an ADRC 
encoding of-the-error function is used. This approach is desirable as it 
is computationally inexpensive. For the process, a list E, consisting of 
the points e* = I yi - x* I is defined. Defining emin and emax as 
respectively the smallest and largest values from the list, then eDR = 
emax-emin. An encoded error value can then defined as: 

gi = (ei - emin)nl/eDR 

where nl represents the number of quantization levels for 
requantizing ej in a similar manner to the ADRC process described 
above. 

Thus, new lists X and Y are generated by selecting only those 
matches where gi is less than some threshold. If the new lists are 
sufficiently long, these lists may be used to generate a refined least 
squares estimate DR r i s . The threshold for gi and the number of 
matches needed before refining the least squares estimation is 
preferably empirically determined. For example, in one embodiment 
for an process involving 8x8x2 horizontally subsampled blocks and nl 
is 10, only matches corresponding to gi = 0 are used, and the estimate 
is refined only when the new lists contain at least 30 matches. 

In an alternate embodiment, DR estimation can be improved by 
clipping potential DR values and recomputing a DR estimate. In 
particular, in one embodiment, a list D is composed of member di 
which contains the DR value that would cause Xj to equal yi. More 
precisely: 

di = 2Q(yi-MIN)/(0.5 + qi ) 

Improvement is seen by clipping each di- That is, 

di' = max(min(UB,di), LB) 
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where DR c i s is then computed to be the average of d**. The clipped 
method (DRd s ) may be combined with other DR estimates , e.g., DRi se 
in a weighted average to produce a final DR value. For example, the 
weighted average DRest is determined according to the following: 

DRest = wi(DRcis) + w 2 (DRise). 

The weights wi and w 2 are preferably empirically determined by 
examining resultant estimations and images generated therefrom 
from particular weightings. In one embodiment wi = 0.22513 and W2 = 
0.80739. 

The invention has been described in conjunction with the 
preferred embodiment. It is evident that numerous alternatives, 
modifications, variations and uses will be apparent to those skilled in 
the art in light of the foregoing description. 
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CLAIMS 

What is claimed is: 

1. A method for recovering lost /damaged attribute data in a 
bitstream of encoded data comprising attribute data and encoded 
sample data, said method comprising the steps of: 

retrieving decoded neighboring data 445; and 

estimating the lost/damaged attribute data using said encoded 

sample data, said decoded neighboring data, and available attribute 

data 465. 

2. The method as set forth in claim 1, wherein the encoded 
data comprises a video signal. 

3. The method as set forth in claim 1, wherein the encoded 
data comprises an audio signal. 

4. The method as set forth in claim 1, wherein the 

lost/ damaged attribute data comprises a dynamic range value (DR) of 
the samples and available attribute data comprises a minimum value 
(MEM) of said samples, wherein said estimating the lost/damaged 
attribute data comprises the steps of: 

gathering square error data of an estimation function with 
regard to said lost/damaged DR; and 

selecting a DR that minimizes said square error data of said 
estimation function. 
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5. The method as set forth in claim 1, wherein 
lost/damaged attribute data comprises a dynamic range value (DR) of 
the samples and available attribute data comprises a minimum value 
(MIN) of said samples, and wherein estimating comprises selecting 
an estimated value (DR') that minimizes the following function: 

f(DR') = l(y i -g(q i/ MIN / DR*)) 2 

i 

where yi is a neighboring decoded sample value, qi is an 
encoded value of said sample, and g() corresponds to a decoding 
operation. 

6. The method as set forth in claim 1, wherein 
lost/damaged attribute data comprises a dynamic range value (DR) of 
the samples and available attribute data comprises a minimum value 
(MIN) of said samples, and wherein estimating comprises selecting an 
estimated value (DR') that minimizes the following function: 

f(DR') = s[ yj -MIN-^^J 

where yi is a neighboring decoded sample value, qi is an 
encoded value of said sample, and Q corresponds to a quantization 
number of said samples. 

7. The method as set forth in claim 1, wherein 

lost/ damaged attribute data comprises a dynamic range value (DR) of 
the samples and available attribute data comprises a minimum value 
(MIN) of said samples, and wherein estimating comprises selecting an 
estimated value (DR) that minimizes the following function: 
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2 



f(DR') = I 



yi -MIN 



where yi is a neighboring decoded sample value, qi is an 
encoded value of said sample, and Q corresponds to a quantization 



number of said samples. 

8. The method as set forth in claim 1, wherein 
lost/damaged attribute data comprises a dynamic range value (DR) of 
the samples and available attribute data comprises a minimum value 
(MIN) of said samples, and said estimating comprises the following 
equation: 



where Q is the quantization number of said samples, yi is a 
neighboring decoded sample value, and qi is an encoded value of said 



9. The method as set forth in claim 1, wherein 
lost/damaged attribute data comprises a dynamic range value (DR) of 
the samples and available attribute data comprises a minimum value 
(MIN) of said samples, and said estimating comprises the following 
equation: 



DR' = 



(2 Q -l)-Z(y i -MIN)-q i 



sample. 



DR' = 



2 Q £( yi -MIN)-q, 
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where Q is the quantization number of said samples, yi is a 
neighboring decoded sample value, and qi is an encoded value of said 
sample. 

10. The method as set forth in claim 1, wherein 
lost/damaged attribute data comprises a minimum value (MIN) of the 
samples and available attribute data comprises a dynamic range (DR) 
of said samples, said method further estimating said the lost/damaged 
attribute data comprising the steps of: 

gathering square error data of an estimation function with 
regard to said lost/ damaged MIN; and 

selecting an estimated MIN value (MIN*) that minimizes said 
square error data of said estimation function. 

11. The method as set forth in claim 1, wherein 
lost/ damaged attribute data comprises said MIN of the samples, 
available attribute comprises said DR of said samples and wherein 
estimating comprises selecting an estimated MIN value (MIN 1 ) that 
minimizes the following function: 

f(MIN> I( yi -g( qi ,MIN',DR)) 2 

i 

where y\ is a neighboring decoded sample value, qi is an 
encoded value of said sample, and g() corresponds to a decoding 
operation. 

12. The method as set forth in claim 1, wherein 
lost/damaged attribute data comprises said MIN of the samples, 
available attribute comprises said DR of said samples and wherein 
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estimating comprises selecting an estimated MIN value (MINT) that 
minimizes the following function: 



where yi is a neighboring decoded sample value, qi is an 
encoded value of said sample, and Q corresponds to a quantization 
number of said samples. 

13. The method as set forth in claim 1, wherein 
lost/ damaged attribute data comprises said MIN of the samples, 
available attribute comprises said DR of said samples and wherein said 
estimating comprises selecting an estimated MIN value (MIN') that 
minimizes the following function: 



where yi is a neighboring decoded sample value, qi is an 
encoded value of said sample, and Q corresponds to a quantization 
number of said samples. 

14. The method as set forth in claim 1, wherein clipping 
comprises limiting estimated attribute data to the range of [Lq, Uq], 
where Lq and Uq respectively represent the lower and upper bounds 
of the range of possible attribute values for Q quantization bits. 

.15. The method as set forth in claim 1, wherein 
lost/ damaged attribute data comprises said MIN of the samples, 
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available attribute comprises said DR of said samples and wherein said 
estimating comprises the following equation: 



MIN' = - 

N 



where Q is the quantization bit of said samples, yi is a 
neighboring decoded sample value, qi is an encoded value of said 
sample, and N corresponds to the number of terms used in the 
summation. 



16. The method as set forth in claim 1, wherein 
lost/ damaged attribute data comprises said MIN as the samples, 
available attribute comprises said DR of said samples and wherein said 
comprises the following equation: 



MIN' = 

N 



where Q is the quantization bit of said samples, y\ is a 
neighboring decoded sample value, qj is an encoded value of said 
sample, and N corresponds to the number of terms used in the 
summation. 



17. The method as set forth in claim 1, wherein 
lost /damaged attribute data comprises a dynamic range value (DR) 
and a minimum value (MIN) of the samples, the method further 
estimating said the lost/damaged attribute data comprising the steps 
of: 

gathering square error data of an estimation function with 
regard to said lost/ damaged DR and said lost/ damaged MIN; and 
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selecting a DR and a MIN that minimize said square error data 
of said estimation function. 

18. The method as set forth in claim 1, wherein 

lost/ damaged attribute data comprises a dynamic range value (DR) 
and a minimum value (MIN) of the samples, and wherein estimating 
comprises selecting estimated values of DR' and MIN' that minimize 
the following function: 

f(DR' ,MIN*) = Z(y. -gfa^MIN' ,DR')) 2 

i 

where yj is a neighboring decoded sample value, qi is an 
encoded value of said sample, and g() corresponds to a decoding 
operation. 

19. The method as set forth in claim 1, wherein 
lost/damaged attribute data comprises a dynamic range value (DR) 
and a minimum value (MIN) of the samples, and wherein estimating 
comprises selecting estimated values of DR' and MIN' that minimize 
the following function: 

f (DR 1 ,MW) = z(y, " MIN 1 J 

where yi is a neighboring decoded sample value, q i is an 
encoded value of said sample, and Q corresponds to a quantization 
number of said samples. 

20. The method as set forth in claim 1, wherein 

lost/ damaged attribute data comprises a dynamic range value (DR) 



Printed from Mimosa 



WO 99/21090 



PCT/US98/22411 



67 



and a minimum value (MIN) of the samples, and wherein estimating 
comprises selecting estimated values of DR' and MIN' that minimize 
the following function: 



f(DR' ,MIN') = - MIN' - qi 2 Q R j 



where yi is a neighboring decoded sample value, qi is an 
encoded value of said sample, and Q corresponds to a quantization 
number of said samples. 



21. The method as set forth in claim 1, wherein 
lost/damaged attribute data comprises a dynamic range value (DR) 
and a minimum value (MIN) of the samples, and wherein estimating 
comprises selecting estimated values of DR' and MIN' that minimize 
the following functions: 



DR = 



(2 Q -1) 



N-Zqi-yj-IqiSyi 



N 



DR r ^ 



MIN'=- 

N 

where Q is the quantization bit number of said samples, yi is a 
neighboring decoded sample value, qi is an encoded value of said 
sample, and N corresponds to the number of terms used in the 
summation. 



22. The method as set forth in claim 1, wherein 
lost /damaged attribute data comprises a dynamic range value (DR) 
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and a minimum value (MIN) of the samples, and wherein estimating 
comprises selecting estimated values of DR' and MIN' that minimize 
the following functions: 



where Q is the quantization bit number of said samples, yi is a 
neighboring decoded sample value, qi is an encoded value of said 
sample, and N corresponds to the number of terms used in the 
summation. 

23. The method as set forth in claim 1, wherein 
lost/damaged attribute data comprises a minimum value (MIN) of the 
samples and available attribute data comprises a dynamic range value 
(DR) of said samples, and wherein said method further comprises 
estimating the lost/damaged attribute data by. calculating a plurality of 
estimated data of the lost/damaged attribute data using said encoded 
sample data, said neighboring decoded data, and available attribute 
data; and selecting the median of said plurality of estimated data. 

24. The method as set forth in claim 1, wherein 
lost/damaged attribute data comprises a minimum value (MIN) of the 
samples, available attribute comprises a dynamic range value (DR) of 
said samples and said estimating comprises selecting the median of a 
plurality of estimated data (MIN') according to the following equation: 





N 
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MIN* = med(MIN , i ) = medj y. 




where yi is a neighboring decoding sample value, qi is an 



encoded value of said sample, Q corresponds to a quantization 
number of said samples, and med() corresponds to median function. 

25. The method as set forth in claim 21, wherein 
lost/damaged attribute data comprises a minimum value (MIN) of the 
samples, available attribute comprises a dynamic range value (DR) of 
said samples and said estimating comprises selecting the median of a 
plurality of estimated data (MIN 7 ) according to the following equation: 



where yi is a neighboring decoded sample value, q i is an 
encoded value of said sample, Q corresponds to a quantization 
number of said samples, and med() corresponds to a median function. 

26. The method as set forth in claim 1, wherein estimating 
said lost/ damaged attribute data using said encoded sample data, said 
neighboring decoded data, and available attribute data, further 
comprises the steps of: 

selecting said decoded neighboring for said estimation according 
to said motion; and 

estimating lost/damaged attribute data using said encoded 
sample data, said selected neighboring decoded data, and available 
attribute data. 
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27. The method as set forth in claim 1, wherein 
lost/damaged attribute data comprises a dynamic range value (DR) of 
the samples and available attribute data comprises a minimum value 
(MIN) of said samples, and wherein estimating said the lost/damaged 
attribute data (DR') comprises the steps of: 

estimating said lost/damaged attribute data using said encoded 
sample data, said selected neighboring decoded data, and available 
attribute data; and 

clipping said estimated lost/damaged attribute data (DR') 
according to the following equation: 

MIN + DR' < NUM . 

28. The method as set forth in claim 27, wherein NUM is 
equal to 255. 

29. The method as set forth in claim 1, wherein 
lost/damaged attribute data comprises a minimum value (MIN) of the 
samples and available attribute data comprises a dynamic range value 
(DR) of said samples, and wherein said estimating said lost/damaged 
attribute data (MIN') comprises the steps of: 

estimating said lost/ damaged attribute data using said encoded 
sample data, said selected neighboring decoded data, and available 
attribute data; and 

clipping said estimated lost/damaged attribute data (MIN') 
according to the following equation: 

MIN' +DR < NUM . 

30. The method as set forth in claim 29, wherein NUM is " 
equal to 255. 



Printed from Mimosa 



WO 99/21090 



PCT/US98/22411 



71 



31. The method as set forth in claim 1, wherein 
lost/ damaged attribute data comprises a dynamic range value (DR) 
and a minimum value (MIN) of the samples, and wherein said 
estimating said lost/damaged attribute data (DR' and MIN') further 
comprises the steps of: 

estimating lost/damaged attribute data using said encoded 
sample data, said selected neighboring decoded data, and available 
attribute data; and 

clipping said estimated lost/damaged attribute data (DR' and 
MIN') according to the following equation: 

MIN'+DR'<NUM. 



32. The method as set forth in claim 31, wherein NUM is 
equal to 255. 

33. The method as set forth in claim 1, wherein said 
estimating said lost/ damaged attribute data using said encoded sample 
data, said neighboring decoded data, and available attribute data, 
further comprises the steps of: 

detecting correlation between said encoded data and 
corresponding neighboring decoded data; 

selecting said neighboring decoded data for said estimation 
according to said correlation; and 

estimating said lost/damaged attribute data using said encoded 
sample data, said selected neighboring decoded data, and available 
attribute data. 
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34. The method as set forth in claim 33, wherein 
lost/ damaged attribute data comprising a dynamic range value (DR) of 
the samples, available attribute data comprising a minimum value 
(MIN) of said samples, is said encoded sample value, yi is said 
neighboring decoded sample value, and wherein detecting correlation 
comprises further steps of: 

defining the range of {Lqi Uqi } corresponding to an encoded 
sample data q\, wherein Lqi is the lower bound corresponding to said 
qi and Uqi is the upper bound corresponding to said qi ; 

comparing a decoded neighboring data yj corresponding to said 
encoded data qi with said range {Lqi Uqi} ; and 

deciding an introduction of said decoded neighboring data yj to 
said estimation step according to said comparison. 



35. The method as set forth in claim 34, wherein 
lost/ damaged attribute data comprising a dynamic range value (DR) of 
the samples, available attribute data comprising a minimum value 
(MIN) of said samples, and said comparison comprising further 
equations: 

(yi-MINjm 
max^ -0.5,0) ~ Q 

( 7i -MIN)m 
min( qi -0.5,0) Q 

where yi is a neighboring decoded sample value, qi is an 
encoded sample value, Q is the number of quantization bits for said 
samples, m represents the maximum quantization code (2Q-1), Lq 
corresponds to the lower bound of Q bits, Uq corresponds to the upper 
bound of Q bits, and max and min respectively corresponds maximum 
and minimum functions. 
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36. The method as set forth in claim 14, wherein the lower 
and upper bounds of the clipping range of attribute values for Q 
quantization bits are relaxed outside the range of possible attribute 
values for suppressing visual degradation caused by instability of 
recovered attribute values. 



37. The method as set forth in claim 36, wherein Lq and Uq 
are respectively adjusted to Lq and U'q according to the following 
equations: 

L'q = r - Lq-i + (1-r) • Lq 
Uq = (1-r) • Uq + r - U Q+1 

wherein r is a relaxation constant, Lq - 1 is the lower bound of 
the range of possible attribute values for Q - 1 quantization bits, and 
Uq+i is the upper bound of the range of possible attribute values for Q 
+ 1 quantization bits. 



38. The method as set forth in claim 37, wherein r is equal to 

.5. 



39. The method as set forth in claim 1, wherein the encoded 
sample data and neighboring decoded data is grouped according to 
relative direction and wherein said estimating comprises: 

generating directional estimates of the lost/damaged attribute 
data for each direction using corresponding encoded sample data and 
encoded neighboring data; 

weighting the directional estimates according relative amounts 
of correlation; and 
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combining the weighted directional estimates to generate 
a combined estimate. 



40. A system comprising a processor configured to recover 
lost/damaged attribute data in a bitstream of encoded data comprising 
attribute data and encoded sample data, said processor configured to: 
retrieve decoded neighboring data 445; and 
estimate the lost/damaged attribute data using said encoded 
sample data, said decoded neighboring data, and available attribute 
data 465. 



41. The system as set forth in claim 40, wherein the encoded 
data comprises a video signal. 



42. The system as set forth in claim 40, wherein the encoded 
data comprises an audio signal. 



43. The system as set forth in claim 40, wherein 
lost/ damaged attribute data comprises a dynamic range value (DR) of 
the samples and available attribute data comprises a minimum value 
(MIN) of said samples, and wherein said processor is configured to 
estimate the lost/damaged attribute data by gathering square error data 
of an estimation function with regard to said lost/damaged DR and 
selecting a DR that minimizes said square error data of said estimation 
function. 



44. The system as set forth in claim 40, wherein 
lost/ damaged attribute data comprises a dynamic range value (DR) of 
the samples and available attribute data comprises a minimum value 
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(MIN) of said samples, and wherein said processor is configured to 
estimate by selecting an estimated value (DR') that minimizes the 
following function: 

f(DR') = I(y, - g^MD^DR')) 2 

where y 4 is a neighboring decoded sample value, qi is an 
encoded value of said sample, and g() corresponds to a decoding 
operation. 



45. The system as set forth in claim 40, wherein 
lost/damaged attribute data comprises a dynamic range value (DR) of 
the samples and available attribute data comprises a minimum value 
(MIN) of said samples, and said processor is configured to estimate by 
selecting an estimated value (DR') that minimizes the following 
function: 

f(I»')-z(y,-MlN-fL^J 

where yj is a neighboring decoded sample value, qi is an 
encoded value of said sample, and Q corresponds to a quantization 
number of said samples. 



46. The system as set forth in claim 40, wherein 
lost/damaged attribute data comprises a dynamic range value (DR) of 
the samples and available attribute data comprises a minimum value 
(MIN) of said samples, and said processor is configured to estimate by 
selecting an estimated value (DR') that minimizes the following 
function: 
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where yi is a neighboring decoded sample value, qi is an 
encoded value of said sample, and Q corresponds to a quantization 
number of said samples. 



47. The system as set forth in claim 40, wherein 
lost/damaged attribute data comprises a dynamic range value (DR) of 

l e T ,!n mpleS ^ 3VaiIable a " ribUte data com Prises a minimum value 
(MIN) of said samples, said processor is configured to estimate using 
the following equation: 



DR'= (2Q ~ l) '? (yi ~ MIN)q ' 

i 

where Q is the quantization number of said samples, yi is a 
neighboring decoded sample value, and qj is an encoded value of said 
sample. 



48. The system as set forth in claim 40, wherein 
lost/damaged attribute data comprises a dynamic range value (DR) of 

^T 163 aV3ilable 3ttribUte d3ta com P™<* a minimum value 
(MIN) of said samples, said processor is configured to estimate using 
the following equation: 

Sq? 
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where Q is the quantization number of said samples, yi is a 
neighboring decoded sample value, and q, is an encoded value of said 
sample. 



49. The system as set forth in claim 40, wherein 
lost/damaged attribute data comprises a minimum value (MIN) of the 
samples and available attribute data comprises a dynamic range (DR) 
of said samples, said system further configured to estimate the 
lost/damaged attribute data by gathering square error data of an 
estimation function with regard to said lost/damaged MIN, and 
selecting an estimated MIN value (MIN') that minimizes said square 
error data of said estimation function. 



50. The system as set forth in claim 40, wherein 
lost/damaged attribute data comprises said MIN of the samples, 
available attribute comprises said DR of said samples and said ' 
processor is configured to estimate by selecting an estimated MIN 
value (MIN') that minimizes the following function: 

f(MIN') = I(y. -g^MnN^DR)) 2 

where yi is a neighboring decoded sample value, qi is an 
encoded value of said sample, and g() corresponds to a decoding 
operation. 



51. The system as set forth in claim 40, wherein 
lost/damaged attribute data comprises said MIN of the samples, 
available attribute comprises said DR of said samples and said ' 
processor is configured to select an estimated MIN value (MIN') that 
minimizes the following function: 
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f(MIN') = Z(y, - MIN* J 



where yj is a neighboring decoded sample value, q, is an 
encoded value of said sample, and Q corresponds to a quantization 
number of said samples. 



52. The system as set forth in claim 40, wherein 
lost/ damaged attribute data comprises said MIN of the samples, 
available attribute comprises said DR of said samples and said 
processor is configured to estimate by selecting an estimated MIN 
value (MIN 1 ) that minimizes the following function: 



f(MIN') = ijy, - MIN J 



where y\ is a neighboring decoded sample value, qi is an 
encoded value of said sample, and Q corresponds to a corresponding 
to a quantization number of said samples. 



53. The system as set forth in claim 40, wherein the processor 
is configured to clip by limiting estimated attribute data to the range of 
[Lq, Uq], where Lq and Uq respectively represent the lower and upper 
bounds of the range of possible attribute values for Q quantization bits. 



54. The system as set forth in claim 40, wherein 
lost/ damaged attribute data comprises said MIN of the samples, 
available attribute comprises said DR of said samples and the 
processor is configured to estimate using the following equation: 
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MIN=-^ 2 -1 ) 

N 

where Q is the quantization bit of said samples. yi is a 
neighboring decoded sample value, qj is an encoded value of said 
sample, and N corresponds to the number of terms used in the 
summation. 



55. The system as set forth in claim 40, wherein 
lost/damaged attribute data comprises said MIN as the samples 
available attribute comprises said DR of said samples and the 
processor is configured to estimate using the following equation: 

MIN' = -^ 2 > 

N 

where Q is the quantization bit of said samples, Tl is a 
neighboring decoded sample value, qi is an encoded value of said 
sample, and N corresponds to the number of terms used in the 
summation. 



56. The system as set forth in claim 40, wherein a 
lost/damaged attribute data comprises a dynamic range value (DR) 
and a minimum value (MIN) of the samples, the processor further 
configured to estimate the lost/damaged attribute data by gathering 
square error data of an estimation function with regard to said 
lost/damaged DR and said lost/damaged MIN, and selecting a DR and 
a MIN that minimize said square error data of said estimation 
function. 
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57. The system as set forth in claim 40, wherein 
lost/damaged attribute data comprises a dynamic range value (DR) 
and a minimum value (MIN) of the samples, and the processor is 
configured to estimate by selecting estimated values of DR' and MIN' 
that minimize the following function: 

f(DR\MIN>i;( yi -g(q„MIN',DR')) 2 

where yi is a neighboring decoded sample value, qi is an 
encoded value of said sample, and g() corresponds to a decoding 
operation. 



58. The system as set forth in claim 40, wherein 
lost/damaged attribute data comprises a dynamic range value (DR) 
and a minimum value (MIN) of the samples and the processor is 
configured to select estimated values of DR' and MIN' that minimize 
the following function: 



f(DR*,MiN')=5:f y^-MiN--^; 01 *'! 2 

' V 2 Q - 1 ) 



where yi is a neighboring decoded sample value, qi is an 
encoded value of said sample, and Q corresponds to a quantization 
number of said samples. 



59. The system as set forth in claim 40, wherein 
lost/damaged attribute data comprises a dynamic range value (DR) 
and a minimum value (MIN) of the samples and the processor is 
configured to select estimated values of DR' and MIN' that minimize 
the following function: 
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f(DR" ,MIN') = l(y, - MIN' J 



where yi is a neighboring decoded sample value, qj is an 
encoded value of said sample, and Q corresponds to a quantization 
number of said samples. 



60. The system as set forth in claim 40, wherein 
lost/damaged attribute data comprises a dynamic range value (DR) 
and a minimum value (MIN) of the samples and the processor is 
configured to select estimated values of DR' and MIN' that minimize 
the following functions: 



MIN 

N 



where Q is the quantization bit number of said samples, yi is a 
neighboring decoded sample value, qi is an encoded value of said 
sample, and N corresponds to the number of terms used in the 
summation. 



61. The system as set forth in claim 40, wherein 
lost/damaged attribute data comprises a dynamic range value (DR) 
and a minimum value (MIN) of the samples and the processor is 
configured to select estimated values of DR' and MIN' that minimize 
the following functions: 
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MIN' = 



DR' = 




N 



where Q is the quantization bit number of said samples, yi is a 
neighboring decoded sample value, qi is an encoded value of said 
sample, and N corresponds to the number of terms used in the 
summation. 



62. The system as set forth in claim 40, wherein 
lost/damaged attribute data comprises a minimum value (MIN) of the 
samples and available attribute data comprises a dynamic range value 
(DR) of said samples, and wherein said processor is further configured 
to estimate the lost/damaged attribute data by calculating a plurality of 
estimated data of the lost/damaged attribute data using said encoded 
sample data, said neighboring decoded data, and available attribute 
data, and selecting the median of said plurality of estimated data. 

63. The system as set forth in claim 40, wherein 
lost/damaged attribute data comprises a minimum value (MIN) of the 
samples, available attribute comprises a dynamic range value (DR) of 
said samples and said processor is configured to estimate by selecting 
the median of a plurality of estimated data (MIN') according to the 
following equation: 
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where y, is a neighboring decoding sample value, qi is an 
encoded value of said sample, Q corresponds to a quantization 
number of said samples, and med() corresponds to median function. 

64. The system as set forth in claim 40, wherein 
lost/damaged attribute data comprises a minimum value (MIN) of the 
amples, available attribute comprises a dynamic range value Sr) of 
axd samples and the processor is configured to estimate by sel ctl 
he med ia n of a plurality of estimated data (MIN') according to the 
following equation: 5 

MIN'= medCMnsJ',.) = med(y., q ' 2 f R j 

where yi is a neighboring decoded sample value, qi is an 
encoded value of said sample, Q corresponds to a quantization 
number of said samples, and med() corresponds to a median function. 

65. The system as set forth in claim 40, wherein the processor 
is configured to estimate said lost/damaged attribute data using said 
encoded sample data, said neighboring decoded data, and available 
attribute data, said processor configured to select said decoded 
neighboring for said estimation according to said motion, and 
estimate lost/damaged attribute data using said encoded sample data 
sa:d selected neighboring decoded data, and available attribute data. ' 

66. The system as set forth in claim 40, wherein 
tost/damaged attribute data comprises a dynamic range value (DR) of 
hesamples .and available attribute data comprises a minimum value 
(MIN) of said samples, the processor is configured to estimate the ' 
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lost/damaged attribute data (DR') by estimating said lost/damaged 
attribute data using said encoded sample data, said selected 
neighboring decoded data, and available attribute data, and clipping 
said estimated lost/damaged attribute data (DR') according to the 
following equation: 

MIN + DR* < NUM . 



67. The system as set forth in claim 66, wherein NUM is 
equal to 255. 

68. The system as set forth in claim 40, wherein 

lost/ damaged attribute data comprises a minimum value (MIN) of the 
samples and available attribute data comprises a dynamic range value 
(DR) of said samples, the processor is configured to estimate the 
lost/damaged attribute data (MIN') by estimating said lost/damaged 
attribute data using said encoded sample data, said selected 
neighboring decoded data, and available attribute data, and clipping 
said estimated lost/damaged attribute data (MIN') according to the 
following equation: 

MIN* +DR < NUM . 



69. The system as set forth in claim 68, wherein NUM is 
equal to 255. 

70. The system as set forth in claim 40, wherein 

lost/ damaged attribute data comprises a dynamic range value (DR) 
and a minimum value (MIN) of the samples, and the processor is 
configured to estimate lost/damaged attribute data (DR' and MIN') by 
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estimating lost/ damaged attribute data using said encoded sample 
data, said selected neighboring decoded data, and available attribute 
data and clipping said estimated lost/damaged attribute data (DR' and 
MIN') according to the following equation: 

MIN*+DR*<NUM. 



71. The system as set forth in claim 70, wherein NUM is 
equal to 255. 

72. The system as set forth in claim 40, wherein the processor 
is configured to estimate said lost/damaged attribute data using said 
encoded sample data, said neighboring decoded data, and available 
attribute data, said processor further configured to: 

detect correlation between said encoded data and corresponding 
neighboring decoded data; 

select said neighboring decoded data for said estimation 
according to said correlation; and 

estimate said lost/damaged attribute data using said encoded 
sample data, said selected neighboring decoded data, and available 
attribute data. 



73. The system as set forth in claim 72, wherein 
lost/ damaged attribute data comprising a dynamic range value (DR) of 
the samples, available attribute data comprising a minimum value 
(MIN) of said samples, qi is said encoded sample value, yi is said 
neighboring decoded sample value, and the processor is configured to 
detect correlation by 

defining the range of {Lqi Uqi} corresponding to an encoded 
sample data qi, wherein Lqi is the lower bound corresponding to said 
qi and Uq* is the upper bound corresponding to said qi; 
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comparing a decoded neighboring data yi corresponding to said 
encoded data qi with said range {Lqi Uqi}; and 

deciding an introduction of said decoded neighboring data y* to 
said estimation step according to said comparison. 

74. The system as set forth in claim 73, wherein 
lost/damaged attribute data comprising a dynamic range value (DR) of 
the samples, available attribute data comprising a minimum value 
(MIN) of said samples, and said comparison comprising further 
equations: 

( 7i -MIN)m ^ L 
max(qj -0.5,0) Q 

(Yt-MIN)m 
min^ -0.5,0) Q 

where yi is a neighboring decoded sample value, qi is an 
encoded sample value, Q is the number of quantization bits for said 
samples, m represents the maximum quantization code (2Q-1), Lq 
corresponds to the lower bound of Q bits, Uq corresponds to the upper 
bound of Q bits, and max and min respectively corresponds maximum 
and minimum functions. 

75. The system as set forth in claim 53, wherein the lower 
and upper bounds of the clipping range of attribute values for Q 
quantization bits are relaxed outside the range of possible attribute 
values for suppressing visual degradation caused by instability of 
recovered attribute values. 
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76. The system as set forth in claim 75, wherein Lq and Uq 
are respectively adjusted to L'q and Uq according to the following 
equations: 

L'q = r • Lq-i + (1-r) • Lq 
Uq = (1-r) • Uq + r . U Q+a 

wherein r is a relaxation constant, Lq . i is the lower bound of 
the range of possible attribute values for Q - 1 quantization bits, and 
Uq+i is the upper bound of the range of possible attribute values for Q 
+ 1 quantization bits. 

77. The system as set forth in claim 76, wherein r is equal to 

.5. 

78. The system as set forth in claim 40, wherein the encoded 
sample data and neighboring decoded data is grouped according to 
relative direction and the processor is configured to estimate by 
generating directional estimates of the lost/ damaged attribute data for 
each direction using corresponding encoded sample data and encoded 
neighboring data, weighting the directional estimates according 
relative amounts of correlation, and combining the weighted 
directional estimates to generate a combined estimate. 

79. A computer readable medium containing executable 
instructions which, when executed in a processing system, cause the 
system to recover lost/damaged attribute data in a bitstream of 
encoded data comprising attribute data and encoded sample data, the 
instructions causing the system to perform the following steps 
comprising : 

retrieving decoded neighboring data 445; and 
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estimating the lost/ damaged attribute data using said encoded 
sample data, said decoded neighboring data, and available attribute 
data 465. 



80. The computer readable medium as set forth in claim 79, 
wherein the encoded data comprises a video signal. 

81. The computer readable medium as set forth in claim 79, 
wherein the encoded data comprises an audio signal. 

82. The computer readable medium as set forth in claim 79, 
wherein lost/ damaged attribute data comprises a dynamic range value 
(DR) of the samples and available attribute data comprises a 
minimum value (MIN) of said samples, and wherein said estimating 
the lost/damaged attribute data comprises executable instructions, 
which when executed in the processing system cause the system to 
perform the following steps: 

gathering square error data of an estimation function with 
regard to said lost/damaged DR; and 

selecting a DR that minimizes said square error data of said 
estimation function. 

83. The computer readable medium as set forth in claim 79, 
wherein lost /damaged attribute data comprises a dynamic range value 
(DR) of the samples and available attribute data comprises a 
minimum value (MIN) of said samples, and said estimating 
comprising executable instructions, which when executed in the 
processing system cause the system to select an estimated value (DR*) 
that minimizes the following function: 
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f(DR') = S(yi - g(q i ,MIN / DR')) 2 

where yi is a neighboring decoded sample value, qi is an 
encoded value of said sample, and g() corresponds to a decoding 
operation. 

84. The computer readable medium as set forth in claim 79, 
wherein lost/ damaged attribute data comprises a dynamic range value 
(DR) of the samples and available attribute data comprises a 
minimum value (MIN) of said samples, and said estimating 
comprising executable instructions, which when executed in the 
processing system cause the system to select an estimated value (DR 1 ) 
that minimizes the following function: 

f(DR') = z[ yi -MIN-^^] 2 

where yi is a neighboring decoded sample value, qi is an 
encoded value of said sample, and Q corresponds to a quantization 
number of said samples. 

85. The computer readable medium as set forth in claim 79, 
wherein lost /damaged attribute data comprises a dynamic range value 
(DR) of the samples and available attribute data comprises a 
minimum value (MIN) of said samples, and said estimating 
comprising executable instructions, which when executed in the 
processing system cause the system to select an estimated value (DR') 
that minimizes the following function: 

f(DR>l(y,-MIN-ad^y 
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where yi is a neighboring decoded sample value, qi is an 
encoded value of said sample, and Q corresponds to a quantization 
number of said samples. 



86. The computer readable medium as set forth in claim 79, 
wherein lost/damaged attribute data comprises a dynamic range value 
(DR) of the samples and available attribute data comprises a 
minimum value (MINT) of said samples, and said estimating 
comprising executable instructions, which when executed in the 
processing system cause the system to perform said estimating using 
the following equation: 

(2 Q -l)-l(y i -MIN).q 1 
DR' = i 5 

i 

where Q is the quantization number of said samples, yi is a 
neighboring decoded sample value, and qj is an encoded value of said 
sample. 



87. The computer readable medium as set forth in claim 79, 
wherein lost/damaged attribute data comprises a dynamic range value 
(DR) of the samples and available attribute data comprises a 
minimum value (MIN) of said samples, and said estimating 
comprising executable instructions, which when executed in the 
processing system cause the system to perform said estimating using 
the following equation: 

2 Q -l(y 1 -MIN)-q i 
DR* = «- 



i 
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where Q is the quantization number of said samples, yi is a 
neighboring decoded sample value, and qi is an encoded value of said 
sample. 

88. The computer readable medium as set forth in claim 79, 
wherein lost/damaged attribute data comprises a minimum value 
(MIN) of the samples and available attribute data comprises a dynamic 
range (DR) of said samples, and said estimating comprising executable 
instructions, which when executed in the processing system cause the 
system to perform the following steps: 

gathering square error data of an estimation function with 
regard to said lost/damaged MIN; and 

selecting an estimated MIN value (MIN') that minimizes said 
square error data of said estimation function. 

89. The computer readable medium as set forth in claim 79, 
wherein lost/damaged attribute data comprises said MIN of the 
samples, available attribute comprises said DR of said samples and 
said estimating comprising executable instructions, which when 
executed in the processing system cause the system to perform a step 
of selecting an estimated MIN value (MIN 1 ) that minimizes the 
following function: 

f(MIN , ) = X(y i -g(q i ,MIN , / DR)) 2 

where yi is a neighboring decoded sample value, qi is an 
encoded value of said sample, and g() corresponds to a decoding 
operation. 
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90. The computer readable medium as set forth in claim 79, 
wherein lost/damaged attribute data comprises said MIN of the 
samples, available attribute comprises said DR of said samples and 
said estimating comprising executable instructions, which when 
executed in the processing system cause the system select an estimated 
MIN value (MIN*) that minimizes the following function: 

f(MINO = l( yi -MIN'-|^J 

where yi is a neighboring decoded sample value, qi is an 
encoded value of said sample, and Q corresponds to a quantization 
number of said samples. 



91. The computer readable medium as set forth in claim 79, 
wherein lost/damaged attribute data comprises said MIN of the 
samples, available attribute comprises said DR of said samples and 
said estimating comprising executable instructions, which when 
executed in the processing system cause the system to perform the step 
selecting an estimated MIN value (MIN*) that minimizes the 
following function: 



f(MIN') = l(y, - MIN' -Si^J 



where yj is a neighboring decoded sample value, qj is an 
encoded value of said sample, and Q corresponds to a corresponding 
to a quantization number of said samples. 



92. The computer readable medium as set forth in claim 79, 
wherein clipping comprises executable instructions, which when 
executed in the processing system cause the system to limit estimated 
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attribute data to the range of [Lq, Uq], where Lq and Uq respectively 
represent the lower and upper bounds of the range of possible 
attribute values for Q quantization bits. 



93. The computer readable medium as set forth in claim 79, 
wherein lost/ damaged attribute data comprises said MIN of the 
samples, available attribute comprises said DR of said samples and 
said estimating comprising executable instructions, which when 
executed in the processing system cause the system to perform the step 
using the following equation: 

mq\t= — — I 

N 

where Q is the quantization bit of said samples, yi is a 
neighboring decoded sample value, qi is an encoded value of said 
sample, and N corresponds to the number of terms used in the 
summation. 



94. The computer readable medium as set forth in claim 79, 
wherein lost/damaged attribute data comprises said MIN as the 
samples, available attribute comprises said DR of said samples and 
said estimating comprising executable instructions, which when 
executed in the processing system cause the system to perform the said 
estimating using the following equation: 

MIN' = -^ ± — — I 

N 

where Q is the quantization bit of said samples, y\ is a 
neighboring decoded sample value, qi is an encoded value of said 
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sample, and N corresponds to the number of terms used in the 
summation. 

95. The computer readable medium as set forth in claim 79, 
wherein a lost /damaged attribute data comprises a dynamic range 
value (DR) and a minimum value (MIN) of the samples, further 
comprising executable instructions, which when executed in the 
processing system cause the system the steps of: 

gathering square error data of an estimation function with 
regard to said lost/damaged DR and said lost/damaged MIN; and 

selecting a DR and a MIN that minimize said square error data 
of said estimation function. 

96. The computer readable medium as set forth in claim 79, 
wherein lost/ damaged attribute data comprises a dynamic range value 
(DR) and a minimum value (MIN) of the samples, 

said estimating comprising executable instructions, which 
when executed in the processing system, cause the system to select 
estimated values of DR' and MIN' that minimize the following 
function: 

f(DR* ,MIN') = - g( qi ,MIN' ,DR')) 2 

where yi is a neighboring decoded sample value, qi is an 
encoded value of said sample, and g() corresponds to a decoding 
operation. 

97. The computer readable medium as set forth in claim 79, 
wherein lost/ damaged attribute data comprises a dynamic range value 
(DR) and a minimum value (MIN) of the samples, 
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said estimating comprising executable instructions, which 
when executed in the processing system, cause the system to select 
estimated values of DR' and MIN' that minimize the following 
function: 



where yi is a neighboring decoded sample value, qiis an 
encoded value of said sample, and Q corresponds to a quantization 
number of said samples. 

98. The computer readable medium as set forth in claim 79, 
wherein lost/damaged attribute data comprises a dynamic range value 
(DR) and a minimum value (MIN) of the samples, 

said estimating comprising executable instructions, which 
when executed in the processing system, cause the system to select 
estimated values of DR' and MIN' that minimize the following 
function: 



encoded value of said sample, and Q corresponds to a quantization 
number of said samples. 

99. The computer readable medium as set forth in claim 79, 
wherein lost/damaged attribute data comprises a dynamic range value 
(DR) and a minimum value (MIN) of the samples, 

said estimating comprising executable instructions, which 
when executed in the processing system, cause the system select 





where yi is a neighboring decoded sample value, qi is an 
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estimated values of DR' and MIN' that minimize the following 
functions: 



DR 

N 



(2 Q -l){N.Iq ryi -I q ,.Z yi ] 



MIN' = - 

N 

where Q is the quantization bit number of said samples, yi is a 
neighboring decoded sample value, qi is an encoded value of said 
sample, and M corresponds to the number of terms used in the 
summation. 



100. The computer readable medium as set forth in claim 79, 
wherein lost/damaged attribute data comprises a dynamic range value 
(DR) and a minimum value (MIN) of the samples, 

said estimating comprising executable instructions, which 
when executed in the processing system, cause the system select 
estimated values of DR' and MIN' that minimize the following 
functions: 



2 Q . 



DR* = - 



MIN 



N-Iq^Iq/] 



N 
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where Q is the quantization bit number of said samples, yi is a 
neighboring decoded sample value, qj is an encoded value of said 
sample, and N corresponds to the number of terms used in the 
summation. 



101. The computer readable medium as set forth in claim 79, 
wherein lost/damaged attribute data comprises a minimum value 
(MEM) of the samples and available attribute data comprises a dynamic 
range value (DR) of said samples, said medium comprising executable 
instructions, which when executed in the processing system, cause the 
system to estimate the lost/damaged attribute data comprising the 
steps of: 

calculating a plurality of estimated data of the lost/damaged 
attribute data using said encoded sample data, said neighboring 
decoded data, and available attribute data; and selecting the median of 
said plurality of estimated data. 

102. The computer readable medium as set forth in claim 79, 
wherein lost/damaged attribute data comprises a minimum value 
(MIN) of the samples, available attribute comprises a dynamic range 
value (DR) of said samples and said estimating comprising executable 
instructions, which when executed in the processing system, cause the 
system select the median of a plurality of estimated data (MEN') 
according to the following equation: 

MIN' = med(MIN , i ) = med^y,- ^ Q ' P * j 

where yi is a neighboring decoding sample value, qi is an 
encoded value of said sample, Q corresponds to a quantization 
number of said samples, and med() corresponds to median function. 



Printed from Mimosa 



WO 99/21090 



PCT/US98/22411 



98 

103. The computer readable medium as set forth in claim 79, 
wherein lost/damaged attribute data comprises a minimum value 
(MIN) of the samples, available attribute comprises a dynamic range 
value (DR) of said samples and said estimating comprises selecting the 
median of a plurality of estimated data (MIN') according to the 
following equation: 

MIN'= med(MIN' i ) = med^ - 

where y\ is a neighboring decoded sample value, qi is an 
encoded value of said sample, Q corresponds to a quantization 
number of said samples, and med() corresponds to a median function. 

104. The computer readable medium as set forth in claim 79, 
wherein estimating said lost/damaged attribute data using said 
encoded sample data, said neighboring decoded data, and available 
attribute data comprises executable instructions, which when executed 
in the processing system, cause the system to perform the steps of: 

selecting said decoded neighboring for said estimation according 
to said motion; and 

estimating lost/ damaged attribute data using said encoded 
sample data, said selected neighboring decoded data, and available 
attribute data. 

105. The computer readable medium as set forth in claim 79, 
wherein lost/ damaged attribute data comprises a dynamic range value 
(DR) of the samples and available attribute data comprises a 
minimum value (MIN) of said samples, estimating said the 
lost/damaged attribute data (DR') comprising executable instructions, 
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which when executed in the processing system, cause the system to 
perform the steps of: 

estimating said lost/damaged attribute data using said encoded 
sample data, said selected neighboring decoded data, and available 
attribute data; and 

clipping said estimated lost/damaged attribute data (DR') 
according to the following equation: 

MIN + DR'< NUM. 



106. The computer readable medium as set forth in claim 105, 
wherein NUM is equal to 255. 

107. The computer readable medium as set forth in claim 79, 
wherein lost/ damaged attribute data comprises a minimum value 
(MIN) of the samples and available attribute data comprises a dynamic 
range value (DR) of said samples, estimating said lost/damaged 
attribute data (MIN') comprising executable instructions, which when 
executed in the processing system, cause the system to perform the 
steps of: 

estimating said lost/ damaged attribute data using said encoded 
sample data, said selected neighboring decoded data, and available 
attribute data; and 

clipping said estimated lost /damaged attribute data (MIN') 
according to the following equation: 

MIN'+DR <NUM. 



108. The computer readable medium as set forth in claim 107, 
wherein NUM is equal to 255. 
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109. The computer readable medium as set forth in claim 79, 
wherein lost/damaged attribute data comprises a dynamic range value 
(DR) and a minimum value (MIN) of the samples, estimating said 
lost/ damaged attribute data (DR' and MUM') further comprising 
executable instructions, which when executed in the processing 
system, cause the system to perform the steps of: 

estimating lost/damaged attribute data using said encoded 
sample data, said selected neighboring decoded data, and available 
attribute data; and 

clipping said estimated lost/damaged attribute data (DR' and 
MIN') according to the following equation: 

MIN'+DR'< NUM. 



110. The computer readable medium as set forth in claim 109, 
wherein NUM is equal to 255. 

111. The computer readable medium as set forth in claim 79, 
wherein estimating said lost/damaged attribute data using said 
encoded sample data, said neighboring decoded data, and available 
attribute data, further comprising executable instructions, which when 
executed in the processing system, cause the system to perform the 
steps of: 

detecting correlation between said encoded data and 
corresponding neighboring decoded data; 

selecting said neighboring decoded data for said estimation 
according to said correlation; and 

estimating said lost /damaged attribute data using said encoded 
sample data, said selected neighboring decoded data, and available 
attribute data. 
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112. The computer readable medium as set forth in claim 111, 
wherein lost/damaged attribute data comprising a dynamic range 
value (DR) of the samples, available attribute data comprising a 
minimum value (MIN) of said samples, qi is said encoded sample 
value, yi is said neighboring decoded sample value, and said detecting 
correlation comprising executable instructions, which when executed 
in the processing system, cause the system to further perform the steps 
of: 

defining the range of (Lqi Uqi} corresponding to an encoded 
sample data qi, wherein Lqi is the lower bound corresponding to said 
qiand Uqi is the upper bound corresponding to said qi; 

comparing a decoded neighboring data yi corresponding to said 
encoded data qi with said range (Lqi Uqi); and 

deciding an introduction of said decoded neighboring data yi to 
said estimation step according to said comparison. 



113. The computer readable medium as set forth in claim 112, 
wherein lost/damaged attribute data comprising a dynamic range 
value (DR) of the samples, available attribute data comprising a 
minimum value (MIN) of said samples, and said instruction 
comprising the comparison comprising further equations: 

(y t -MIN)m 
max(q i -0.5,0) Q 

( yi -MIN)m ^ u 
min(q. -0.5,0) Q 

where yi is a neighboring decoded sample value, qi is an 
encoded sample value, Q is the number of quantization bits for said 
samples, m represents the maximum quantization code (2Q-1), Lq 
corresponds to the lower bound of Q bits, Uq corresponds to the upper 
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bound of Q bits, and max and min respectively corresponds maximum 
and minimum functions. 



114. The computer readable medium as set forth in claim 92, 
wherein the lower and upper bounds of the clipping range of attribute 
values for Q quantization bits are relaxed outside the range of possible 
attribute values for suppressing visual degradation caused by 
instability of recovered attribute values. 

115. The computer readable medium as set forth in claim 114, 
wherein Lq and Uq are respectively adjusted to Lq and U'q according 
to the following equations: 

Lq= r • Lq-i + (1-r) • Lq 
U'q = (1-r) - Uq + r • U Q+ i 

wherein r is a relaxation constant, Lq . i is the lower bound of 
the range of possible attribute values for Q - 1 quantization bits, and 
Uq+i is the upper bound of the range of possible attribute values for Q 
+ 1 quantization bits. 

116. The computer readable medium as set forth in claim 115, 
wherein r is equal to .5. 

117. The computer readable medium as set forth in claim 79, 
wherein the encoded sample data and neighboring decoded data is 
grouped according to relative direction and wherein estimating 
comprises executable instructions, which when- executed in the 
processing system, cause the system to perform the steps of : 
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generating directional estimates of the lost/damaged attribute 
data for each direction using corresponding encoded sample data and 
encoded neighboring data; 

weighting the directional estimates according relative amounts 
of correlation; and 

combining the weighted directional estimates to generate 
a combined estimate. 



118. An apparatus configured to recover lost/ damaged 
attribute data in a bitstream of encoded data comprising attribute data 
and encoded sample data, said apparatus comprising: 

means for retrieving decoded neighboring data; and 

means for estimating the lost/damaged attribute data using said 

encoded sample data, said decoded neighboring data, and available 

attribute data. 

119. The apparatus as set forth in claim 118, wherein 

lost/ damaged attribute data comprises a dynamic range value (DR) of 
the samples and available attribute data comprises a minimum value 
(MIN) of said samples, wherein said apparatus comprising means to 
estimate the lost/damaged attribute data by gathering square error data 
of an estimation function with regard to said lost/damaged DR and 
selecting a DR that minimizes said square error data of said estimation 
function. 



120. The apparatus as set forth in claim 118, wherein 
lost/damaged attribute data comprises a minimum value (MIN) of the 
samples and available attribute data comprises a dynamic range (DR) 
of said samples, wherein said apparatus comprising means for 
estimating the lost/damaged attribute data by gathering square error . 
data of an estimation function with regard to said lost/damaged MIN, 
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and selecting an estimated MIN value (MEN') that minimizes said 
square error data of said estimation function. 

121. The apparatus as set forth in claim 118, wherein said 
apparatus estimates said lost/damaged attribute data using said 
encoded sample data, said neighboring decoded data, and available 
attribute data, said apparatus further comprising: 

means for detecting correlation between said encoded data and 
corresponding neighboring decoded data; 

means for selecting said neighboring decoded data for said 
estimation according to said correlation; and 

means for estimating said lost/ damaged attribute data using 
said encoded sample data, said selected neighboring decoded data, and 
available attribute data. 



122. The apparatus as set forth in claim 118, wherein the 
encoded sample data and neighboring decoded data is grouped 
according to relative direction and the apparatus comprises means to 
estimate by generating directional estimates of the lost/damaged 
attribute data for each direction using corresponding encoded sample 
data and encoded neighboring data, means for weighting the 
directional estimates according relative amounts of correlation, and 
means for combining the weighted directional estimates to generate a 
combined estimate. 
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MEASUREMENTS: SQUARE ERROR 



S. E-ZEfo-yjj)' 
1 J 



where 



Xj is the (candidate) decoded value of pixel i 
yy is a horizontal or vertical neighbor of Xj 



*NOTE: Current implementation discards 3 largest terms, i.e., (Xj - yjj) 2 



FIG. 20a 
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